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16511 GROUND MOVEMENTS CAUSED BRACED EXCAVATIONS 


KEY WORDS: Construction procedure; Deformation; Excavation; 
Foundations; Lateral pressure; Loads (forces); Performance evaluation; 
Settlement control; Soil-structure interaction 


ABSTRACT: The causes and patterns ground movement associated with braced 
excavations are examined. Field measurements displacements for braced cuts 
sand and soft medium clay show relationship between the deflection the 
excavation wall and the ratio horizontal vertical movement adjoining the cut. 
review various case histories indicates that soil displacement depends strongly 
the construction process. Practices are recommended for promoting stiff bracing 


system, excavating beneath the lowest braces, and using berms restrain ground 
movement. 


REFERENCE: O’Rourke, Thomas (Asst. Prof., School Civ. and Environmental 
Engrg., Cornell Univ., Ithaca, N.Y. 14853), “Ground Movements Caused Braced 
Excavations,” Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. 
GT9, Proc. Paper 16511, September, 1981, pp. 1159-1178 


16510 LATERALLY LOADED PILE DESIGN 


KEY WORDS: Analysis (mathematics); Coefficient subgrade reaction; 
Computer analysis; Computers; Earth pressure; Foundations; Lateral forces; 
Piles; Piles (foundations); Soil mechanics; Soil-structure interaction 


ABSTRACT: beam-bar finite element formulation laterally loaded pile problem 
developed. The approach taken utilizes the beam elastic foundation coefficient 
subgrade reaction method. The three loading conditions lateral load applied 
moment free-end pile and lateral load applied fixed-end pile are analyzed 
for range soil strength variations with depth. Curves are developed which present 
the maximum values function the soil-structure stiffness ratio. design 
procedure with example given show the use these curves the design 
any laterally loaded pile system. 


REFERENCE: Sogge, Robert (Principal Consulting Geotechnical Engr., Desert 
Earth Engrg., Tuscon, Ariz. 85705), “Laterally Loaded Pile Design,” Journal the 
Geotechnical Engineering Division, ASCE, Vol. 107, No. GT9, Proc. Paper 16510, 
September, 1981, pp. 1179-1199 


16493 SEISMIC RESPONSE SOFT CLAY SITES 


KEY WORDS: Clays; Computer programs; Cyclic loads; 
Damping; Degradation; Earthquakes; Mathematical models; San Francisco; 
Seismic stability; Shear modulus; Site evaluation; Soft soils; Vibration 
response 


ABSTRACT: nonlinear stress-strain model, utilizing Ramberg-Osgood idealization 
backbone curve and allowing for its progressive degradation due cyclic loading, 
used conduct site response studies two soft clay sites. The sites selected are: (1) 
Southern Pacific Building, San Francisco, and (2) Station No. 119 Site, Tokyo, Japan. 
The subsurface conditions and the material properties used analyses are presented. 
The computer program DCHARM used for nonlinear analysis allows for progressive 
degradation the backbone curve. For comparison the results obtained from program 
SHAKE using equivalent linear idealization are also presented. The program 
DCHARM shown provide ground response results closer the recorded data 
than SHAKE. The results obtained indicate that substantially higher accelerations and 
durations earthquakes are needed cause significant degradation soft clay. 


REFERENCE: Singh, Ram (Sr. Project Engr., Woodward-Clyde Consultants, 5120 
Butler Pike, Plymouth Meeting, Pa. 19462), Debry, Ricardo, Doyle, Earl H., and 
Idriss, Izzat M., “Nonlinear Seismic Response Soft Clay Sites,” Journal the 
Geotechnical Engineering Division, ASCE, Vol. 107, No. GT9, Proc. Paper 16493, 
September, 1981, pp. 1201-1218 


' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 
' 


ce 


16491 ANISOTROPIC STRENGTH COHESIONLESS SANDS 


KEY WORDS: Anisotropic soils; Anisotropy; Cohesionless soils; Fabrics; 
Granular materials; Sands; Sandy soils; Shear strength; Soil mechanics; Soil 
strength 


ABSTRACT: Natural soil deposits are generally anisotropic well heterogeneous. 
order provide guideline for estimating reasonable strength parameter 
applicable such complex soils, plane strain tests uniform-anisotropic, uniform- 
pseudoisotropic, layered-anisotropic, and layered-pseudoisotropic samples have been 
performed. Anisotropic strength response, caused the anisotropic distribution 
contact normals, should expected even uniform mass spheres. always 
necessary assume that sand generally anisotropic. Tests samples prepared 
the plunging-tapping method are advisable for obtaining reasonable estimate the 
mean strength for anisotropic sands. 


REFERENCE: Oda, Masanobu (Visiting Scholar, Dept. Civ. Engrg., Northwestern 
Univ., Evanston, Ill.; (Per, Address) Dept. Foundation Engrg., Saitama Univ., 
Japan), “Anisotropic Strength Cohesionless Sands,” Journal the Geotechnical 
Engineering Division, ASCE, Vol. 107, No. GT9, Proc. Paper 16491, September, 
1981, pp. 1219-1231 


16489 GEOTEXTILE-REINFORCED UNPAVED ROAD DESIGN 


KEY WORDS: Aggregates; Fabrics; Geotextiles; Reinforced roads; 
Reinforcement; Road design; Soil stability; Surface treatment (roads); 
Textiles; Unpaved roads 


ABSTRACT: providing reinforcement, geotextiles improve the performance 
unpaved roads, and for given traffic, the thickness the aggregate layer can 
reduced. This paper explores the mechanism geotextile reinforcement roadways, 
and method developed which enables the engineer calculate the required 
thickness aggregate layer and make proper selection the geotextile used. 
Charts are presented assist the design geotextile-reinforced unpaved roads and 
evaluate the reduction aggregate thickness resulting from the use geotextile. 
Practical examples and comparisons with full scale tests are presented. 


REFERENCE: Giroud, Jean-Pierre(eWoodward-Clyde Consultants, Geotextiles and 
Geomembranes Group, Chicago, 60603), and Noiray, Laure, 
Reinforced Unpaved Road Design,” Journal the Geotechnical Engineering Division, 
ASCE, Vol. 107, No. GT9, Proc. Paper 16489, September, 1981, pp. 1233-1254 
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GROUND MOVEMENTS CAUSED 
BRACED EXCAVATIONS 


INTRODUCTION 


The design braced excavations performed the basis two fundamental 
requirements: (1) The selection excavation and support measures necessary 
for stability the cut; and (2) the control ground movements affecting 
adjacent property. Stability generally means preventing failure, i.e., designing 
the braced cut remain open without danger construction personnel 
very large movements that jeopardize surrounding structures. Design measures 
for stability include the selection earth pressures, determining the potential 
for base soil failure, and the layout dewatering system prevent piping 
and blowout along the walls and bottom the cut. The control ground 
movements generally based the types buildings and public facilities 
near the excavation. Foremost among the measures necessary for control 
the prediction displacements various distances from the cut. addition, 
control measures include the protection adjacent structures and the use 
construction methods restrict ground movements levels consistent with 
the acceptable performance surrounding facilities. 

With regard stability, apparent earth pressure envelopes are recommended 
for variety soil types (16,19,20) and methods are available for estimating 
the likelihood plastic deformation and base failure for excavations clay 
(5,16). contrast, the patterns settlement associated with braced cuts are 
summarized only very general terms, and horizontal displacements are treated 
mostly case history basis. The wide range performance resulting from 
different soil conditions and construction methods often complicates the deci- 
sion-making process regarding both the potential for damage and protection 
nearby structures. 

This paper examines the ground movements caused braced excavations 
with emphasis the relationship between ground movements and various aspects 
the construction process. Sources lost ground are reviewed and deformation 
patterns the excavation wall are related displacement patterns the ground 

Prof., School Civ. and Environmental Engrg., Cornell Univ., Ithaca, N.Y. 
14853. 

Note.—Discussion open until February 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication October 31, 
1980. This paper part the Journal the Geotechnical Engineering Division, Proceedings 


the American Society Civil Engineers, ©ASCE, Vol. 107, No. GT9, September, 
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surface adjoining the excavation. Case histories braced cuts both sands 
and soft clays are summarized and recommendations for ground movement 
control are made the basis observational data. 


Site 


Braced cuts commonly are visualized step-wise process excavation 
and support. Following each increment excavation, inward deflections 
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Construction sequence 


FIG. Displacement Record for (14 Deep Cut Soft Medium 
Clay 


the walls occur. These deflections cause movements the soil adjoining the 
cut that are transmitted the ground surface settlements and horizontal 
offsets. This process forms the conceptual basis for much the analysis and 
finite element modeling braced cut performance and has been helpful for 
understanding the ground response actual excavation and bracing procedures. 
Nevertheless, deep excavations often are performed combination with other 


activities that cause movement particularly during the preparatory stages 
construction. 


ft 
(43 
(30) 
(61) 
= | | Install concrete piles \ 
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Site preparation may include: (1) Relocation and underpinning utilities; 
(2) dewatering aquifers above and below the base excavation; (3) construction 
the excavation wall; and (4) the installation deep foundations. some 
cases, the movements associated with site preparation will exceed those occurring 
result the excavation and support process. The relocation utilities, 
e.g., may have locally severe impact adjacent property, especially when 
trenching carried close other pipelines and communication conduits. 
Alternatively, dewatering may consolidate the soil over area that exceeds 
substantially the area affected excavation-induced movements. Wall construc- 
tion may require predrilling for soldier piles, the use vibratory hammers 
install sheet piles, the excavation slurry panels for concrete diaphragm 
walls. Each these can cause permanent movements, the magnitude and 
distribution which will vary according the soil conditions and details 
the construction procedures. 

Ground movements associated with site preparation are shown the records 
deep basement construction. example, Fig. relates typical settlement 
and horizontal movements with the various construction activities performed 
for (14 deep braced cut the Embarcadero III Project San 
Francisco. The cut was excavated through approximately (6.5 granular 
fill and sands and soft medium silts and clays. The bottom 
the cut was underlain roughly (11 medium clay followed 
deposit interbedded dense sand and stiff clay. The general construction 
procedure this site described Drossel (7). The soil conditions and temporary 
bracing system were similar those associated with braced cut construction 
the Levi Strauss Building, which considered Tait and Taylor (18). 

The ground movements were strongly influenced construction activities 
related site preparation for deep excavation and pile driving. Over 800 timber 
piles were removed from the site and similar number prestressed concrete 
piles were driven. Sheet piles were installed along the perimeter the cut 
with vibratory hammer. Between Day and Day 170, approximately 30% 
the total ground movements occurred even though the excavation was 
maintained constant (3.4 depth. some cases, there was clear 
relationship between movement and construction activity, e.g., when timber 
piles were removed near the property line, horizontal movements were induced 
that cracked the street surface for distance ft-25 (4.5 m-7.6 
from the excavation’s edge. 

similar manner, the construction drilled shafts Chicago has led 
large settlements. Fig. summarizes the settlements adjacent three braced 
cuts where large movements were caused, principally the result ground 
loss during drilled shaft, caisson, construction. The cases are selected 
show extreme conditions ground loss and are not necessarily typical general 
construction practices. The settlements and distances are plotted dimensionless 
form fractions the maximum excavation depth. Zones settlement, 
delineated Peck (18) for various soil types and excavation conditions, are 
shown for comparison. 

Chicago practice for drilled shaft construction has been considered detail 
Baker and Khan (1) and ground movements associated with the construction 
method have been described Lukas and Baker (10). Therefore, only brief 
description given here. general, oversized shafts were drilled depths 
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ranging from ft-35 m-10.5 below the construction surface. After 
temporary steel linings were set, the shafts were further extended augering 
and advancing steel casing, which was left place part the permanent 
foundation. All drilled shafts were carried the limestone bedrock, approximately 
110 (34 below street level. 

Much the settlement associated with case resulted from dewatering 
difficulties deposit sand, silt, and gravel directly overlying bedrock. 
Pumping was required prepare the drilled shaft bottoms for concrete and, 
the process, sand and silt were pumped out with the water. Large portions 
the settlements for Case and Case were caused clay squeezing into 
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Maximum depth excavation 
1.0 2.0 


after Peck (16) 
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FIG. 2.—Settlements Adjacent Braced Cuts with Caisson Construction Difficulties 


= 
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the shafts during excavation and into annular voids around the shafts after 
steel linings had been installed. all cases, the movements associated with 


drilled shaft construction were between 50% and 70% the total settlements 
shown Fig. 


Braced excavations generally are performed three prominent stages. During 
each stage, excavation and support methods contribute the ground displacement 
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characteristic way. These stages include: (1) Initial excavation before bracing; 
(2) excavation subgrade after the upper braces are installed; and (3) removal 
the braces. 

Fig. traces the development lateral soil movement function the 
construction history for cut-and-cover excavation Washington, D.C. The 
excavation was (18 deep and was supported soldier pile and lagging 
walls with five levels cross-lot struts. The soldier piles 130) were 
(23 long and were located 7.5 (2.3 centers. Field observations 
this excavation and details the instrumentation have been reported 
the writer (12). 

The figure shows the horizontal wall movements and cumulative horizontal 
strains the soil for each the construction stages. The horizontal strains 
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FIG. 3.—Horizontal Strains Associated with Various Stages Braced Cut Construc- 
tion 


have been estimated from inclinometer measurements dividing the differential 
lateral displacement between two points given elevation the distance 
separating them and plotting the strain the midpoint between the two 
measurements. addition, vectors cumulative ground movement are shown 
various locations. Each vector was determined from the vertical and horizontal 
movement common point. The horizontal movements were measured with 
inclinometers, and the vertical movements were measured either optical 
leveling the street surface vertical extensometers depth. The contour 
interval for horizontal strain 0.5 which equivalent the tensile 
strain which cracking induced masonry structures (2). Strain magnitudes 
prefixed with minus sign indicate tension. The construction sequence and 
lateral movements are described under the following headings. 


40 
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60 
(183) 
80 
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j 
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Initial Excavation before Bracing.—The excavation was deepened 
m-7.5 before braces were installed. Thus, deformation the wall 
occurred primarily cantilever-type movement. The horizontal strains reflect 
this mode deformation triangular pattern contours that decrease 
magnitude with depth and distance from the wall. 

Excavation Subgrade after Upper Braces braces were installed, 
the upper portion the wall was restrained from further lateral movement. 
fact, preloading the upper level struts decreased the horizontal movement 
the top the cut, resulting slight recompression the soil. the 
deeper portion the cut, inward bulging the wall caused tensile strains, 
the contours for which were inclined approximately 45° from the vertical. 

Removal the bottom braces were removed build the under- 
ground structure, further inward bulging the wall occurred. When the upper 
braces were removed, the wall was supported its lower portion the subway 
Structure resulting cantilever-type deformation the top the wall. 
Consequently, the cumulative strains are composite the horizontal distortion 
associated with inward bulging depth and the distortion associated with 
cantilever movement the upper levels excavation. 

important recognize the plastic deformation the base the cut 
will restrained during brace removal owing the weight the invert 
basement slab that constructed before brace removal begins. For cuts 
soft medium clay, therefore, this construction stage likely result 
only small fraction the ground movements affected base heave during 
initial excavation and excavation subgrade. 

summary, the lateral movements generated open cutting are related 
directly the mode deformation the excavation wall which, turn, 
related the construction procedure. The final movement pattern composed 
essentially two zones horizontal strain. deep-seated zone horizontal 
strain develops response inward bulging the excavation wall and forms 
the lower boundary the influenced soil mass. upper zone horizontal 
strain develops response either cantilever movement lateral translation. 
The two principal modes wall deformation are shown Fig. 3(d). Because 
the soil strains are closely related the modes wall deformation, would 
useful develop quantitative relationship between the displacement pattern 
the wall and movement the ground surface adjacent the cut. This 
relationship explored the following section. 


The measured displacements associated with several braced cuts were studied 
comparing the deformations the excavation walls with the ratios horizontal 
movement settlement the ground surfaces adjoining the cuts. Fig. 
summarizes the observed movements for two braced cuts Chicago that are 
typical the trend shown all the excavations studied. Inclinometer measure- 
ments wall displacement and optical surveys settlement and lateral 
movements were used develop the comparisons. all cases, the field data 
pertain maximum displacements measured after the excavations had been 
taken subgrade and significant portion the braces removed. The data 
relate movements measured distances 0.35-1.0 times the maximum 
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excavation depths from the cuts. Within these distances, the ratios horizontal 
movement settlement were found relatively constant, showing little 
variation function distance from the excavation. 

All data were screened accordance with two guidelines: (1) Each ratio 
lateral vertical displacement derives from the combined measurement 
lateral movement and settlement the same observation point; and (2) measure- 
ments equal less than 1/4 in. mm) were neglected minimize the 
influence survey error the computed ratio. Furthermore, the computed 


Lateral displacement, in. (mm) 
(50) (107) (152) 
2 4 6 


Brace 
levels 


ie} 
10 15 20 25 30 
Ratio of horizontal to vertical movement, 4/4, 


Depth, ft (m) 


Lateral wal! 
Average = 0.60 movements 
Standard deviation = 0.08 


a) Movements related principally to deep inward bulging 


Lateral displacement, in. (mm) 
(50) (102) (152)(203) 
2 a 6 8 


Frequency 


0 OS 10 15 20 25 30 


Depth, ft (m) 


i ; 40 
Ratio of horizontal to vertical movement, d.,/d, (12-2) 


Average = 1.32 movements 
Standard deviation = 0.50 


b) Movements related principally to cantilever deflection 


FIG. 4.—Ratios Measured Horizontal Movements Settlements Ground Surface 
Adjoining Braced Cuts 


tatios were analyzed statistically assuming both normal and log normal 
distributions the data. The assumed distributions gave similar values the 
mean. The values Fig. are for normal distributions which gave slightly 
higher, more conservative ratios horizontal vertical movement. 

There are marked differences the ratios that depend the pattern 
wall movement. Fig. shows the wall deflection for excavation that was 
supported with stiff braces early stage construction. The wall movement 
occurred principally deep inward bulging. The average value horizontal 
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vertical movement adjoining the cut was 0.60. contrast, Fig. 4(b) shows 
the wall deflection for excavation where substantial cantilever movements 
occurred owing lack early support and braces without adequate stiffness. 
The average value horizontal vertical movement adjoining the cut was 

illustrate the influence wall distortion soil movements, coefficient 
the deformation, defined Fig. The numerator the expression 
measure the cantilever movement, the wall. The denominator 
represents the combined cantilever, and inward bulging, components 
displacement. must emphasized that the coefficient used only 
approximate measure relative deformation. intended show general 
trends the data pertaining the immediate displacements caused open 
cutting. Furthermore, the coefficient defined for conditions which the 
bottom the excavation wall secured relatively firm stratum. Consequently, 
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Deformation for Fixed-End Walls 


> 
~ 
x 
v 
2 
ra 
3 
a 
€ 
} 
= 
© 
© 
° 
a 


= 
> 
° 
2 
Cc 
= 
° 
° 
= 
° 
a 


the data summarized Fig. pertain excavations where the walls were 
driven stiff soil layer, and thus, were not affected significant movement 
the wall bottoms. 

Fig. shows the ratio horizontal displacement settlement the ground 
surface adjoining the cut function the coefficient deformation. Average 
values the ratio are plotted. Field observations analyzed terms the 
coefficient were chosen for cases where there was little influence movements 
from either site preparation consolidation. The table the figure summarizes 
information about the braced cuts. Data were derived from total seven 
braced cuts, ranging depth from ft-60 m-18 m). 

The field data show limits for the ratio horizontal vertical movement 
equal approximately 0.6 and 1.6 for wall deformation caused solely inward 
bulging and solely cantilever movement, respectively. These bounds compare 
favorably with the limits determined from model tests, performed Milligan 
(11), which also are plotted Fig. the model tests, soil strains and 
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displacements medium dense sand were measured response wall 
movements caused incremental excavation. The combined field and model 
test data show that cantilever-type movements are allowed dominate during 
construction, the relative magnitude horizontal ground movement increases 
significantly. Conversely, the wall firmly braced early stage 
excavation, inward deflection the wall will lead horizontal movements 
that are significantly less than settlements. 

The data summarized Fig can used estimate horizontal movements, 
either directly from settlement surveys planning basis prior construction. 
Most excavation and bracing schemes will generate both cantilever displacement 
and deep inward movement the wall. Consequently, value one will 
provide reasonable, preliminary estimate the ratio horizontal movement 
settlement. Settlements can estimated the basis local experience 
from field measurements summarized the literature. 

Because settlements and horizontal movements are linked directly with the 
wall movement, important understand how wall deformation influenced 
the excavation and support process. Many factors can distinguished among 
the methods excavation and support that affect ground movements. Three, 
particular, deserve attention. They include brace stiffness, the depth 
excavation beneath the lowest braces, and the use berms. Each considered 
the following sections and examples are taken from field observations 
show ground movement response typical construction procedures. 


Brace 


the level brace contact, wall movements are caused compression 
the interfaces among members the support system, horizontal bending 
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FIG. 6.—Preloading Set-Up and Measured Braced Stiffnesses 


the wale, and elastic the brace. Accordingly, the effective 
stiffness the brace may substantially less than its ideal elastic stiffness. 
Palmer and Kenney (15) report that the effective brace stiffness determined 
from field measurements during open cutting the Oslo Subway was low 
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1/50 the idea elastic stiffness. Studies braced cut behavior Jaworski 
(9) and Scott (17) show that significant horizontal displacements occur when 
braces lack tight connections with the wall compressible materials are used 
shim separations between the support members. 

Commonly, braces are preloaded percentage their design load induce 
rigid contacts among members the support sytem. O’Rourke and Cording 
(13) have considered preloading practices during construction the Washington, 
D.C. Metro. useful examine these practices emphasize some important 
characteristics the preloading operation. 

shown Fig. preloading was performed inserting hydraulic 
jack each side individual pipe strut between the wale and special 
jacking plate welded the strut. Hydraulic pressure was applied until the strut 
was stressed approximately half its design load after which the separation 
between the strut and wale was shimmed with in. mm) thick steel plates. 

Increased separation between the struts and wales were measured during 
preloading. addition, lateral displacements the soldier piles the level 
preloading were measured opposite sides the excavation with inclinome- 
ters attached the soldier piles. The soldier pile displacements were subtracted 
from the measured increases separation obtain the apparent deformation 
each strut. The effective stiffness, strut given level support 


which the average preload; and the average apparent deformation 
each strut level. 


The ideal elastic stiffness, K,, individual strut 


which the modulus elasticity the strut material; and and 
the cross-sectional area and length the strut, respectively. 

determining the ratio effective ideal strut stiffness, possible 
express the brace stiffness fraction its maximum value. Fig. 
shows the ratio effective ideal strut stiffness plotted with respect preload. 
The data were taken from measurements three separate strut levels. The 
relatively high fraction ideal stiffness achieved each strut level shows 
that preloading was effective developing rigid support. 

Although the measurements indicate that strut stiffness increases with increased 
load, very high preloads are not necessarily advantageous. some instances, 
individual struts were preloaded higher than adjacent braces consequence 
variations the preloading technique. the excavation continued, not 
only did the loads these braces remain higher than those adjacent members 
but, some instances, increased faster rate. Very high preloads create 
locally stiff sections the excavation wall that may attract relatively high 
pressures additional excavation occurs. 

Fig. shows soldier pile displacements (shaded area) caused preloading 
level cross-lot struts. The measurements are accurate about 0.01 in. 
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(0.25 mm) the level preloading. The outward movement the soldier 
piles shows that preloading was effective closing separations within the bracing 
system. Movement the level the preloaded struts was maximum 0.1 
in. (2.5 mm) each side the excavation, and was distributed primarily within 
the elevation the jacking and wedging operation. 

Under these conditions, preloads are not likely regain previous soil move- 
ments because they have only localized effect wall displacement. preloads 
are increased above the magnitude needed develop stiff bracing system, 
their contribution the restraint movement becomes increasingly less effective. 
Very high preloads may lead strut loads that exceed design values and may 
cause high bending moments the wall. 


Lateral (mm) 


FIG. 7.—Wall Displacements from Preloading Struts 


There are some instances where substantial movements can occur the brace 
level even though preloading used. Fig. shows relatively common bracing 
procedure which rakers are installed between the wall and portion the 
completed foundation. raker transmitting loads only one two caissons 
may lack the restraint prevent significant wall movement. Elastic analyses, 
such those outlined Davisson and Gill (6), show that the lateral movements 
simple raker and caisson bracing system can high in. (76 mm) 
for braced cut loads typical construction conditions Chicago. Displacements 
this nature lead cantilever deformation the wall. 

Fig. 8(a) shows idealized view wall deformation associated with horizontal 
loading caisson. grade beams are joined several caissons and foundation 
slabs are poured, the restraint transmitted through the raker becomes increasingly 
more efficient. Fig. shows additional wall movement after caisson deflections 
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have ceased and the excavation has been completed. Although the volumes 
incremental displacement shown Figs. 8(a) and are approximately 
equal, the influence wall movement the pattern soil displacement 
substantially different. Fig. the cantilever deformation the wall may 


Horizontal wall 
movement for low 
brace stiffness 


Inward bulging 
wall 


Caisson stabilized 


FIG. 8.—Wall Deflections Associated with Raker and Caisson Bracing 


induce horizontal ground movements that are nearly three times large 
those generated the inward bulging shown Fig. 


Excavation Lowest Braces 


Fig. shows the lateral wall movements for concrete diaphragm wall that 
was used during construction part the Washington, D.C. Metro. The 
(21 wide excavation was extended through approximately (14 
dense sands and stiff clay and (4.5 hard fissured clay. The 


beam 
WA 
Caisson 
Lateral deformation caisson 
Raker 
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water table was located approximately (10.5 below the ground surface. 
The wall was (25 deep and in. (0.85 wide and was reinforced 
with steel piles 124) centers. All braces were preloaded 
approximately 75% their design loads. 

The figure shows lateral wall displacements corresponding three successive 
levels excavation beneath the second level struts. All movements and 
excavation depths are referenced the time when the second level struts were 
installed and preloaded (Day 0). common for cut-and-cover construction, 
the excavation was carried its deepest level the center the cut, leaving 
elevated sections soil against the wall. the central portion the cut 


Lateral displacement, (mm) 


1.2(30) 08(20) 


2nd level Dense sands 
struts and 


stiff clay 


at 
- 
a 


Hard clay 


23, 74- elapsed time 
(days) since installation 
second level struts 


Horizontal Very dense 


(6m) 


FIG. 9.—Observed Lateral Displacements Concrete Diaphragm Wall 


was deepened from ft-55 (10.5 m-17 m), the inward movement increased 
maximum 0.7 in. (18 mm) and the volume displacement nearly doubled. 

The movements shown Fig. emphasize the importance excavation 
depth beneath the lowest braces for controlling ground loss during open cutting. 
The section stiffness the concrete diaphragm wall, which the product 
the wall’s elastic modulus and moment inertia, was approximately 
times that soldier pile and lagging walls used elsewhere during construction 
the Washington, D.C. Metro. Nevertheless, the inward movements the 
diaphragm wall are comparable those reported for the other cuts (12,14). 
The deflections excavation walls under distributed earth loads depends not 
only the séction stiffness, but approximately the fourth power the 


8-sto 
(6) 
(12) 
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(4.5 beneath the lowest braces will offset the benefit derived from 
tenfold increase wall section stiffness. 

When concrete diaphragm walls are used alternative underpinning 
adjacent structures and dewatering, they must support foundation and hydro-static 
pressures that exceed the earth pressures generally imposed soldier pile 
and lagging systems. Under these conditions, unlikely that the deflection 
concrete diaphragm wall will substantially less than that soldier 
pile and lagging wall subject only earth pressure. For good control ground 
movements, experience Washington, D.C. suggests that excavation depths 
limited maximum (5.5 beneath the lowest braces for deep 
cuts medium dense sands and interbedded sands and stiff clay. This limit 
should apply for excavations supported both concrete diaphragm walls and 
soldier pile and lagging systems. 

must emphasized that guidelines this nature not apply for braced 
excavations soft medium clay where conditions base stability influence 
movements beneath the lowest braces. Design measures for conditions large 


plastic deformation the base the cut are considered Peck (16), Clough 
al. (5), and Goldberg al. (8). 


Use 


berm mass soil that left place against the walls excavation. 
Because braced cuts must performed stages, berms lend themselves 
the excavation and support process ways straightforward that they often 
become integral part the construction technique. When deep basements 
are built deposits soft medium clay, berms frequently are used 
restrain the excavation walls the basement structure assembled the 
center the cut. similar manner, berms are used during cut-and-cover 
excavation. the cut extended beneath the lowest braces, the excavation 
may deepened the center remove additional soil and provide path 
for construction vehicles. This practice leaving elevated sections soil adjacent 

Sloughing, seepage, creep, and construction activity will lead the gradual 
distortion berms. Accordingly, instructive look for field examples 
that show how progressive deformations berms have affected braced cut 
performance. 

Fig. shows the inward wall movements for (25 deep braced 
excavation sands and interbedded stiff clay. The braced cut was part 
the construction Gallery Place Station the Washington, D.C. Metro. Details 
the construction and field instrumentation are described elsewhere (12). The 
excavation walls were composed (27.5 long piles 
130) 7.5 (2.3 centers with timber lagging. The struts were prestressed 
half their design loads. Wall displacements and levels excavation are 
referenced the time which the fourth level struts were installed and preloaded 
(Day 0). 

The displacements show the progressive wall movement soil was excavated 
from depth ft-82 (18 m-25 m). The total incremental volume 
movement per unit length given the area between the lateral displacement 
profiles for Day and Day 160. Day 44, approximately half the total incremental 
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movement had occurred response excavation the center the cut. 
From Day Day 147, berm was left place against the wall the 
cut. The approximate dimensions the berm are shown the cut. The 
approximate dimensions the berm are shown the figure. Owing general 
labor strike, additional excavation was performed during this time interval. 
The wall movements increased from Day Day 147, until 80% the total 
incremental displacement had taken place. 
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FIG. 10.—Observed Lateral Displacements and Strut Loads Associated with Berm 
Performance 


Fig. shows the change average load for the fourth level struts 
function time since their installation and preloading. The load increased steadily 
the central portion excavation was deepened with the berm place. 
The average load increased approximately 180 kips (800 kN) from initial 
preload kips (400 kN). Although increased temperature affected the load, 
the influence temperature was small relative that the excavation and 
deformation the berm. 
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The close relationship between deformation and load emphasizes the important 
effects that berms may have the earth pressure distributed the lowest 
braces. Although the magnitude increased displacement caused berm 
deformation was small, the load imposed the bottom struts was very high. 

When berms are used during excavation soft medium clays, the magnitude 
and time dependency movements will increase substantially comparison 
those which will occur for excavation sand. example, Fig. 


summarizes the displacements measured during the excavation basement 
Chicago. 
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FIG. 11.—Observed Displacements for Raker and Berm Excavation Soft Medium 
Clay 


The deep excavation was extended through approximately 
the cut was underlain soft medium clay. The excavation 
wall was composed steel sheet piles (MZ-38) driven stratum stiff 
clay approximately (14 below the ground surface. Measurements are 
summarized for three distinct stages excavation and support. shown 
the figure, these stages include: (1) Excavation the center the cut with 
berm adjacent the wall (2) partial excavation the berm and installation 
the upper rakers; and (3) continued excavation the berm and installation 
the bottom rakers. All displacements and levels excavation are referenced 
the start open cutting (Day 0). 
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shown Fig. 11, the initial excavation was performed conjunction 
with berm, (4.5 wide its top and inclined grade IV:1.5H 
from ft-26 (4.5 m-8 below the ground surface. From Day Day 
26, the volume inward wall movement doubled response gradual 
deformation the berm even though additional excavation was performed. 
Day 45, when the upper level rakers were installed, the berm had been 
cut back about 70% its initial volume. The rakers were not preloaded. 
During this time, wall movements were influenced the driving sheet piles 
for elevator pit. The piles were driven inside the excavation 
for (19.5 distance parallel the north wall. The wall deformation 
was primarily the form cantilever movement, with maximum displacement 
in. (127 mm) the top. the berm was cut back and the bottom rakers 
installed, the wall movements occurred increasingly inward bulging. Fig. 
the maximum wall displacement shown in. (178 mm) the base 
the excavation. The settlements and horizontal movements measured the 
street adjacent the cut are summarized Fig. 11(d). 

The major part the ground movements occurred prior Day 45, when 
the upper rakers were installed. shown increased wall deformation from 
Day Day 36, significant portion these displacements was caused 
gradual yielding the berm adjacent the wall. The measurements 
this site corroborate those Clough and Davidson (3) for excavation 
San Francisco where similar time-dependent loss ground was caused 
berm deformation. 

These observations emphasize the potential for large displacements when berm 
and raker excavations are performed soft medium clays. Successive 
reductions the weight the berm lead increased plastic heave the 
base the excavation. Cutting back the berm, therefore, results loss 
both horizontal and vertical restraint. Under these circumstances, seems prudent 
use large berms and brace the excavation wall firmly early stage 
construction. 

Guidelines for proportioning berms soft medium clay have been developed 
Clough and Denby (4) from finite element simulations. Experience with berm 
and raker construction Chicago indicates that these guidelines can provide 
estimate soil movement function berm size. The time dependency 
berm deformation caused creep and gradual attrition from construction 
activities will contribute soil displacement. Furthermore, the displacement 
associated with berm deformation likely composed relatively large 
component horizontal movement. For these reasons, guidelines and rules 
thumb should used carefully and linked with construction procedures 
that minimize the time lag between excavation the center the cut and 
the installation stiff braces. 


Because braced excavations often call for the integration several distinct 
activities within the general construction plan, ground losses from each activity 
should evaluated when estimating displacements. During site preparation, 
ground movements may caused by: (1) Relocating and underpinning utilities; 
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(2) dewatering; (3) construction the excavation wall; and (4) installation 
deep foundations. 

Field observations and model test data.show direct relationship between 
wall deformation and the ratio horizontal movement settlement adjacent 
the cut. cantilever-type movements are allowed dominate during 
construction, the proportion horizontal vertical movement increases and, 
the limit, approaches value approximately 1.6. Conversely, the wall 
braced firmly early stage excavation, inward deflection the wall 
will lead horizontal movements that are significantly less than settlements. 
Under these conditions, the ratio horizontal movement settlement may 
low 0.6. The ratio will decrease further consolidation occurs. 

Preloading essential promote stiff bracing system and thereby restrict 
wall movement. Nevertheless, practices that call for very high preloads are 
likely cause local concentrations earth pressure that may exceed the pressures 
assumed during design. most cases cross-lot struts prestressed 50% 
their design load will sufficiently rigid restrict further movement the 
level support and sufficiently low load avoid being overstressed 
additional excavation occurs. 

Limiting the excavation depth beneath the lowest braces important for 
preventing ground loss during open cutting. For close control ground move- 
ments experience Washington, D.C. implies that excavation depths should 
limited maximum (5.5 beneath the lowest braces for deep 
cuts medium dense sands and interbedded sands and stiff clay. 

Sloughing, seepage, creep, and construction activity will lead the gradual 
distortion berms. Time-dependent deformations contribute ground loss, 
particulary when berms are used for excavations soft medium clays. 
Experience Chicago indicates that guidelines developed Clough and Denby 
(4) for proportioning berms can provide good estimate soil movement 
function berm size. control ground movements, the time lag between 
excavation the center the cut and the installation stiff braces should 
minimal. 
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The following symbols are used this paper: 
cross-sectional area brace; 


coefficient deformation; 
horizontal ground surface displacement; 
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ground surface settlement; 

elastic modulus brace material; 
effective strut stiffness; 

ideal elastic strut stiffness; 

length brace; 

average preload; 

cantilever wall movement; 

inward bulging wall; and 

average apparent strut deformation. 
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LATERALLY LOADED PILE DESIGN 


INTRODUCTION 


Essentially three approaches (based two different characterizations soil 
strength) have been taken the past analyze laterally loaded piles. One 
the most common methods formulates and solves the differential equation 
for deflection the pile using beam elastic foundation approach (4). 
The soil strength characterized using coefficients subgrade reaction. 
closed form solution the governing differential equation possible when 
the soil modulus constant with depth. Otherwise, nondimensional solutions 
obtained using difference equation formulation must used. 

elasticity approach taken Spillers and Stahl (8) and Poulos (5) describes 
the soil and structure elastic continuum. Such formulation poses increased 
difficulty when the soil modulus varies with depth. The analysis simplified 
neglecting the horizontal shear stresses between the soil and the sides 
the pile. The Poulos solution limited its assumption that the soil modulus 
constant with depth. 

The finite element method can used obtain solutions comparable 
the differential equation and elasticity solutions. the soil strength defined 
using coefficients subgrade reaction, then the difference equation solution 
obtained (2). instead, the soil modeled elastic continuum, the 
elasticity solution approximated (1). 

The method gives 
comparable results solution which idealizes the soil elastic continuum. 
comparison between elastic and subgrade reactions solutions, Poulos (5) 
showed that generally, for all pile flexibilities, the subgrade reaction method 
slightly overestimates the moments, displacements, and rotations. Elastic moduli 
are most often obtained from triaxial compression tests performed the 
laboratory. Coefficients subgrade reaction are determined the field 
horizontal plate loading tests. both, the actual stress path loading should 
simulated. The uncertainties predicting the elastic modulus and the 
coefficient subgrade reaction are similar. 
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this study finite element approach with subgrade reaction soil charac- 
terization used develop curves which eliminate the complex presentations 
which often are associated with other methods. expressing the behavior 
laterally loaded piles common soil-structure interaction approach, clarifica- 
tion the solution fashion useful engineers results. 


Finite Element Model.—In certain mechanics problems the configuration 
the system and its loading such that the soil strength can defined from 
coefficient subgrade reaction relations between load and deformation for 
similar small-scale systems. laterally loaded pile represents such situation. 
Thus coefficients subgrade reaction relating stress displacement rather 
than strain elastic modulus can used. this approach simplified 
one-dimensional bar (spring) soil resistance representation (Fig. can used 
obtain the same results achieved more complex models which represent 
the soil elastic continuum. 


16 beams 
“~ 17 bors 


Note’ For n=O no variation 
with depth occurs 


FIG. 1.—Beam-Bar 
Model 


Finite Element 


FIG. 2.—Net Soil Resistance 


laterally loaded pile can handled using the finite element method 
modeling the pile structure with beam elements and using bars represent 
the uniaxial soil resistance. This approach effect assumes independent 
Winkler-type action for the springs. Any standard beam beam-bar finite element 
program can used design laterally loaded pile systems. bar can 
represented using beam element either setting its moment inertia, 
equal and supporting all three degrees freedom its free end 
setting equal very small number and supporting the translational degrees 
freedom. The equilibrium, stress-strain, and geometric compatibility relations 
for the beam elements are combined form the stiffness relations relating 
the applied loads, the displacements, for each element (6). The stiffness 
relations for bar elements are similarly developed. 

For equilibrium bar element, the internal element forces, are related 
the bar (soil) stresses, 
F=Ap. 


which the contributory contact area the soil. Stress-strain, p-e, relations 
for the bar are 
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which the bar element length. The relation between pressure and 
displacement the system level given 


which the coefficient subgrade reaction. the length the bar 
taken unity Eq. then Eqs. and show that can replaced 
The stiffness relation for bar becomes 


The constitutive equations for this model consist coefficients formed 
the addition the beam element stiffnesses and bar element stiffnesses. For 


0.0 


k=Kmax 


FIG. Reaction Variation with Depth 


these Winkler type spring resistances, the bar element stiffnesses form diagonal 
stiffness array. When added the beam element stiffnesses they form the 


following relation between the applied external nodal loads, and the nodal 
displacements 


[K] {x} 
The example finite element formulation for two-beam and three-bar element 
system presented Appendix 


The finite element model the system presented Fig. The 17-bar 
and 16-beam elements provide adequate division piles any length 
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order that the solution corresponds exactly the closed form differential equation 
solution for constant coefficient subgrade reaction. The model does not 
consider the vertical shear stresses which may exist the pile-soil interface 
nor those the horizontal direction the sides the pile. uniform stress 
distribution across the pile width inherent the method. 

The magnitude and shape distribution the moment along the pile essentially 
identical throughout the pile for the three bottom tip support conditions of: 
(1) Pinned; (2) fixed against rotation and vertical movement; and (3) free except 
for vertical movement. The latter condition has been used this study. For 
cases studied this paper, elastic modulus, and moment inertia, are 
constant along the pile length. The finite element model can also readily 
used develop soil-structure relations for the step-tapered piles often used 
offshore construction, and for those embedded nonhomogeneous soils. 
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FIG. Load Free-End Pile 


Soil Strength Characterization.—The coefficient subgrade reaction represents 
the relation between the net soil pressure, and displacement, for the soil 
surrounding the pile. This net resistance the resisting, passive type stress, 
minus the imposing, active type stress existing opposite sides the pile. 
symmetrical about both the and axis shown Fig. for various 
depths below the surface. The soil stress each side the pile can determined 
multiplying the displacement times coefficients subgrade reaction which 
describe the passive active type contribution the net resistance. These 
values would greater and less, respectively, than the coefficient describing 
the net resistance. The computed changes stress each side the pile 
would additions deletions from the initial at-rest stress each side 
the pile which exists before the lateral load applied. 

The coefficient horizontal subgrade reaction, can related constant 
subgrade reaction, k,, and the pile width, the equation 
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which the vertical distance along the plie measured from the ground 
surface (9). The maximum value occurs when the pile length, giving 


The different variations with depth which can assume are dependent the 
exponent some these are shown Fig. 
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FIG. Load Free-End Pile 


The following values represent secant relation between pressure and 
deflection which useful linear analyses. For soil reactions less than 
ultimate reactions, initial tangent values sufficiently approximate 
secant values. incremental nonlinear analysis would require tangent values 
(Fig. 2). 

Overconsolidated clays, e.g., caused glaciation, exhibit constant value 
values are given 
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300 
kips per cubic foot; 


cubic meters 


which the undrained compressive strength, kips per square foot. 


The lower values are associated with soft clays, the higher values with stiff 
clays. 


0.5 1.0 


FIG. Load Free-End Pile 


Sands usually display linear increase stiffness with depth, i.e., 
The values are the range 


159 


cubic meters 


The maximum value limits the coefficient approximately 100 times 
the passive earth pressure any depth Such limitation valid shallow 
pile depths. deeper pile depths better upper limit would (50 B-100 
kips/cu ft. 

Values the parameters and are dependent the soil type, ground 
water table elevation, drainage conditions, time dependency, and nonlinearity 
present. There large range variation these values since they cannot 
uniquely related these specific conditions. general, submerged soils 
will have values which are less than those for unsubmerged soils. The 
situation dry soil residing over submerged soil can represented 
values which are greater than one. the loading the soils, either 
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due the time rate loading the drainage conditions, approaches 
undrained case, values are closer those given Eq. 11. Soils for 


which the strength rapidly increases with confining pressure have values 
less than one. 


Son-Structure 


Solutions for the three loading conditions free-end pile subject 
lateral load moment and fixed-end pile subject lateral load can best 
presented form common soil-structure interaction problems. this 
form the response parameters, moment, shear, and deflection used describe 
the behavior the pile and the soil pressure are given function the 


0.6 


0.4 


0.2 


0.0 
0.0 1.0 


logS 
FIG. Load Free-End Pile 


ratio, the stiffness the soil that the structure. For pile-soil 
system 


The parameter the ratio soil stiffness structural stiffness and identical 
the nondimensional constant obtained the solution the governing 
differential equation for beam elastic foundation Hetenyi (3). 
this expression since the soil stiffness given for the entire structural section 
should the moment inertia for the section width 

Maximum values moment shear pressure 
and deflection (x,,,,) can plotted for various values Rowe (7) 
presented moment results for anchored sheetpiles similar manner. Limiting 
values for these response parameters can easily calculated for the extreme 
condition rigid pile weak soil. For such case the calculation 


0.3 
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these parameters determinate problem since the deflection pattern known. 
Appendix presents these results. 

For the three loading conditions with piles general flexibilities, interaction 


been developed using the finite element model described previously. These curves 


0.5 1.0 
logS 


FIG. Load Free-End Pile 


0.5 1.0 
logS 


FIG. Load Free-End Pile 


are presented Figs. 4-9 for the loading condition free-end pile subject 
lateral load, Figs. 10-15 for free-end pile subject applied moment, 
and Figs. 16-19 for fixed-end pile subject lateral load. The curves are 
restricted working loads for which the soil strength can characterized 
linear relation. Such working loads are generally the range 1/3-1/2 


105 
10° 
10° 
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the ultimate values. The curves are developed for soil modulus variation 
with depth described shape exponent, equal 0.3, and 

This development was conducted using specific value lateral load 
moment piles, both, idealized Fig. and having different lengths, 


1.0 


FIG. Moment Free-End Pile 


1.0 
logS 


FIG. Moment Free-End Pile 


section properties, and strengths, and embedded soils having different stiff- 
nesses indicated various values coefficients subgrade reaction. 
nondimensional presentation permits simpler development which not all 
the variables need changed. The nondimensional output parameters have 
values which are reduced increased from those for rigid pile. The behavior 
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piles can classified according their relative stiffness ratio into the 
following approximate groupings: 


Rigid: log 0.3. 
Intermediate: 0.3 log 1.3. 
Flexible: log 1.3. 


For relatively rigid piles the solutions given Appendix for perfectly rigid 
piles are applicable. 

Interaction principles are evident these curves (Figs. 4-19). pile can 
made behave relatively rigid manner either being extremely stiff 
itself being embedded extremely weak soil. The more rigid pile 
is, the greater subject larger distances from the top the 


10! 


10° 
0.0 


log 
FIG. Moment Free-End Pile 


pile. For such case, the soil response, measured small and 
occurs deeper depths. 


The following steps should employed the moment design laterally 


loaded pile. The steps this procedure utilize the curves presented Figs. 
4-19. 


Choose trial design pile having specific length, elastic modulus, 
and allowable stress, 

Determine the maximum value the coefficient subgrade reaction, 
and its variation with depth, For overconsolidated clays, may 
zero and given Eqs. and 11. For sands, approximately one 
and given Eq. 12. Any selection these parameters for pure clay 
sand soils well for mixtures each will have some degree uncertainty. 
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FIG. Moment Free-End Pile 
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Develop structural curve for piles various cross section having the 
specific and chosen step This curve determined calculating 
the maximum allowable moment the structure can sustain 


which the distance from the centroid section the extreme outer 
fiber. The moment inertia for reinforced concrete section should based 


FIG. Load Fixed-End Pile 


logS 


FIG. 17.—p,,,,-Lateral Load Fixed-End Pile 


upon the uncracked gross section. Plot the nondimensional parameter incorporat- 
ing this value versus log the appropriate interaction curve. The proper 
design cross section handle the imposed moment given sections above 
the intersection the two curves (Fig. 20). 
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Determine the approximate location for the applicable value 
log 

Determine the values and Evaluate whether 
the design assumption that the maximum net soil pressure, less than 
1/3 1/2 the ultimate unit bearing capacity the soil verified. not, 
the pile width should increased. Approximate ultimate soil pressures are 
q,, for cohesive soils and times the passive Rankine pressure for piles 
cohesionless soils. The maximum soil pressure the passive side the pile 
will greater than the amount the pressure existing the active 
side the pile. 

Determine the approximate locations and p,,,, determined step 
The location and x,,,, always the top the pile. 


1.0 
logS 


FIG. Load Fixed-End Pile 


Design Example.—The design procedure will applied the design 


pile embedded medium dense sand and subject applied k(17.8 
kN) lateral load. 


round timber pile will chosen; 1.5 10° ksi (10.3 MPa), 

For the medium dense sand, (15 20)/B 300/B kef, 


(meters) (0.152) (0.254) (0.356) 

inches* 63.6 201 491 1,020 1,890 
(meters* (0.265) (42.5) (78.7) 

kip-feet 3.53 8.38 16.4 28.3 44.9 
(kilonewton-meters) (4.80) (38.5) (61.1) 


0.0 
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log 


max 


Fig. presents plot superimposing this data Fig. The intersection 
yields 8-in. section, which subject moment kip-feet (11 kN-m). 


Xmax 


FIG. Load Fixed-End Pile 


max 


Timber 
Structural Curve 
(Allowable Mmax) 


Imposed 


0.0 
0.0 0.5 1.0 


logS 
FIG. and Interaction Curves-Design Example 
Fig. M,,,, located 3.5 (1.1 below top pile. 


statics: (17.8 KN); Fig. 2.1 ksf (100 kPa); Fig. 
0.022 radians; and Fig. x,,,, 0.088 (27 mm). 


0.5 
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and x,,,, are located top pile; and Fig. located 
2.5 (0.8 from top pile. less than 1/2 .120 kips/cu 
2.5 2.25 kips/sq (110 kPa). 


For the design example presented, were desired limit 1.5 
ksf (72 kPa), curve’’ giving the maximum soil pressure allowed 
for each pile width versus the log provided that section must developed 
follows: 


(meters) (0.152) (0.203) (0.254) (0.305) (0.356) 

log 1.06 0.940 0.843 0.763 0.696 

3.75 5.00 6.25 7.50 8.75 


structural curve presented Fig. for this design example. The intersection 
the limiting structural curve with the curve, provides lower 


For Pmax 


(86 kPa) 


0.5 
logS 


FIG. 21.—Limiting Design Example 


bound the required pile width. Thus for the design example 10-in. (0.254-m) 
wide pile section greater required. similar procedure could used 
limit this design example. 

the design procedure, structural curve can developed for shear 
similar manner that for moment using instead The structural 
curve for shear, applicable for the case applied moment the top 
the pile, presented Fig. 10. 

Alternatively, design approach which incorporates the pile length, 
variable follows developing set structural curves for piles different 
lengths step the design procedure. For the design problem, using additional 
pile lengths and (3.0 and 12.2 m), the following structural 
curve results: 


Imposed 
0.0 
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(meters) (0.152) (0.254) (0.305) (0.356) 


(3.0 0.088 0.210 0.410 


0.4 


0.5 1.0 


FIG. 22.—Length Selection: Design Example 
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FIG. 23.—Superposition Loadings 


1.12 


0.281 


0.395 
0.997 


Fig. shows such curves. The intersections the structural curves with 
the imposed curve yield the design section which must used for 


each particular pile length. 


For piles having ground surface and top which not coincide, the loading 
can broken into the two cases shown Fig. 23. Since linear has been 
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assumed, the solutions given Figs. 4-9 for lateral load and Figs. 10-15 
for applied moment the ground surface can superimposed. the 
superposition, must remembered that maximums occurring different 
points along the pile cannot added directly. 


Summary 


This study has developed the solution laterally loaded pile formulation 
terms soil-structure interaction procedures. These procedures can 
implemented any beam beam-bar finite element program. 

The finite element method, incorporating linear subgrade reaction soil 
characterization, has been used develop curves defining the maximum moment, 
shear, pressure, deflection, and rotation laterally loaded pile versus the 
soil-structure stiffness ratio. Companion curves which prescribe the location 
these maximum values are developed well. The three tip-end conditions 
of: (1) Lateral load applied free end; (2) moment applied free end; and 
(3) lateral load with fixed-end rotation have been analyzed. Since linear material 
properties have been used, these solutions can superimposed. Four variations 
the coefficient subgrade reaction are evaluated order that bracketing 
soil conditions practical problems possible. should emphasized 
that the solution obtained contains degree uncertainty which due 
predominantly the characterization the soil strength. 

simple design procedure employing the charts outlined. This procedure 
enables the designer quickly obtain the proper design cross section for any 
pile configuration soil having wide range properties. The design procedure 
applied design example where the inherent simplicity the method 
obtaining both the magnitude and location the maximum moment, shear, 
pressure, deflection, and rotation shown. 


formulation developed for the finite element model consisting two-beam 
elements equal length, and three-bar elements shown Fig. 24. The 
stiffness relations for one beam element having the internal forces, and 
deformations, shown Fig. are 


The stiffness relation for individual bar element (Fig. 25) has been given 


F, Sym X4 
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Eq. The internal deformations for the beam and bar elements are identical 
the deformation occurring the nodes the system and thus both are 
denoted 


The system stiffness array for the 5-element model formed adding the 
element stiffnesses representing the internal bar and beam element forces 


BEAM 


FIG. 24.—Finite Element Model For- FIG. Force Deformation 
mulation Notations 


each node. The result the following explicit form Eq. for the 
model: 


+ 


Any given loading condition the ground surface specified assigning 
values the lateral load, applied moment, P,, rotational displacement, 
both, the top the pile. Solution the simultaneous equations for 
the loading condition specified yields the nodal displacements. The internal forces 
the pile can then obtained using Eq. 16, representing the beam element 
stiffness. The soil stress obtained using the nodal deformations Eq. 
for each bar element. 


P, X6 
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For rigid pile embedded soil having stiffness relation 
subjected lateral load magnitude the ground surface, the following 
nondimensional relations exist: 


M max M max 


which obtained from the solution 
n+2 n+l 


(at top) 


Br 


Pmax 


max max 


For this case shifts from near the top the bottom the rigidity becomes 
more pronounced. The value given its magnitude near the top. 

For rigid pile embedded soil having stiffness expressed 
and subject applied moment ground surface the following 
nondimensional relations exist: 


n+2 


1198 


SEPTEMBER 1981 


max 


The following nondimensional relations exist for rigid pile embedded 
soil having stiffness relation with the top the pile being 
fixed the ground surface and subject lateral load, this level: 


max 


max 
(at bottom) 
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The following symbols are used this paper: 


cross-sectional area bar element; 
pile width diameter; 
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elastic modulus; 

internal element force; 

moment inertia; 

array stiffness coefficients; 

coefficient subgrade reaction; 

constant horizontal subgrade reaction; 

element length; 

pile length; 

bending moment; 

shape exponent; 

nodal loads; 

net pressure (stress) soil; 

lateral pile load; 

undrained compressive strength; 

nondimensional soil-structure stiffness ratio 
shear; 

lateral displacement; 

vertical distance along pile measured from ground surface; 
distance from centroid section extreme outer fiber; 

rotation angle radians; and 

stress. 


n 
qu 
x 
Zz 
€ 


Subscripts 


allowable; and 
max maximum. 
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NONLINEAR SEISMIC RESPONSE Sort CLAY SITES 


Members, ASCE and Izzat Idriss,* ASCE 


INTRODUCTION 


The behavior soils under cyclic loading basically nonlinear and hysteretic. 
The soil stress-strain model most widely used analytical procedures for site 
response analysis the equivalent linear technique incorporating strain-dependent 
stiffness and damping. The development this equivalent linear approach has 
included number comparisons between analytical results and recorded ground 
motions. Nonlinear analyses are especially appropriate for soft clay sites where 
nonlinear and modulus degradation effects are more significant than sites 
consisting stiffer soils. However, the use the nonlinear response analysis 
relatively new development, and there need for comparisons between 
motions computed nonlinear methods and calibration with available records. 
This paper summarizes the results such effort. 


There are three basic issues the development realistic analytical procedures 
for site response analyses: (1) Characterization the input motion; (2) modeling 
the soil profile; and (3) stress-strain relation for the soil. Most developments 
and applications assume that the input motions are vertically propagating plane 
shear waves. For horizontally layered sites, response calculation reduces 
the one-dimensional analysis shear beam soil. The computer program 
SHAKE developed Schnabel, al. (22) incorporates the equivalent linear 
model. number case histories have presented comparisons between analytical 
results using this model and recorded ground motions (7,20). Also, some analytical 
studies have been performed comparing results equivalent linear and nonlinear 
calculations (2,8,13,25). 

During the past few years, number nonlinear models stress-strain 
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behavior have been used, including the bilinear (8,17) and multilinear (11) 
representations for the backbone curve. Hyperbolic (18) and Ramberg-Osgood 
(2,9,24) idealizations have also been utilized. number nonlinear site response 
analyses using these models have been presented the literature 
Recently, the computer code DCHARM (9) was developed, which specifically 
designed for the seismic response analysis profiles containing soft clay layers. 
The code DCHARM was obtained modifying the existing program CHARSOIL 
developed Streeter, al. (24), which based the use the method 
characteristics solve the one-dimensional wave equation and uses the 
Ramberg-Osgood model for the backbone curve. The main modification consisted 
adding the capability for degrading the backbone curve the soil during 


cyclic loadings thus incorporating the modulus degradation model for soft clay 
developed Idriss, al. (10). 


For the present study, two sites were selected where strong motion earthquake 
records were available soft clay deposits. These sites are: (1) Southern Pacific 
Building Site, San Francisco, California (March 22, 1957 earthquake); and (2) 
Station No. 119 Site, Tokyo, Japan (July 1968 earthquake). Horizontal 
accelerograms obtained the two sites, having peak accelerations about 
0.05 and total record lengths sec, were used calibrate the 
nonlinear response method allowing for modulus degradation (DCHARM). The 
intensity shaking both cases was high enough for the stress-strain behavior 
the soil become nonlinear, but not high long enough cause significant 
modulus degradation the soft clay. 

the two selected sites, information subsurface conditions and properties 
the soils was available. Nonlinear and degradation parameters for the San 
Francisco Bay Mud, which constitutes the soft layer one the sites, were 
available from earlier cyclic testing program (10). Finally, accelerograms 
recorded simultaneously neighboring rock stiff sites were available both 
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FIG. 1.—Recorded Acceleration Spectra FIG. 2.—Soil Profile and Parameters for 
Rock, San Francisco 
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cases, and were used provide the input excitation. San Francisco, record 
rock was obtained Golden Gate Park about five miles from the soft Southern 
Pacific Building site. Tokyo, record stiff site Station No. 118 
was obtained about miles from the soft site. 


Paciric San Francisco 


earthquake magnitude 5.3 along the San Andreas Fault was recorded 
soil the basement the Southern Pacific Building, San Francisco 
(10.5 miles from the fault) March 22, 1957. The earthquake was simultaneously 
recorded rock the Golden Gate Park USGS (United States Geological 
Survey) station. The Golden Gate Station 7.0 miles from the epicenter and 
5.0 miles from the Southern Pacific site (7). The acceleration response spectra 
the S80E horizontal component recorded Golden Gate for 2.5% and 10.0% 
spectral damping ratios are shown Fig. 


The soil conditions the Southern Pacific building site were obtained from 
previous study (7) and supplemented appropriate values unit weight 
and shear wave velocity for the Franciscan shale rock underlying the site, 
shown Fig. sandy fill about thick covers the site. Below the 
fill 35-ft layer soft clay followed stiff clay down depth 
225 ft; the stiff clay interbedded with 10-ft layer dense sand depth 


125 ft. The stiff clay underlain layer dense sand and gravel, 
which extends bedrock depth 285 ft. The basic cyclic shear stress-strain 
behavior the soil layers was specified the Ramberg-Osgood equation 
the backbone curve the first cycle loading (10): 


Gs Yy Yy 


which: shear stress; shear strain; shear modulus very 
small strains and soil parameters. The experience the 
writers with modulus reduction curves granular soils and stiff very stiff 
clays indicates that, the absence test data, the use the average curve 
for sands proposed Seed and Idriss (19), more appropriate than other 
proposed curves. Therefore, Ramberg-Osgood parameters, 1.0, 
and 0.02%, selected fit the average curve for sands, were used for 
all layers except the 30-ft soft San Francisco Bay Mud layer. For all SHAKE 
runs this study the modulus reduction and damping curves Seed and 
Idriss (19) were used. Cyclic triaxial tests were conducted this soil and 
the results are available elsewhere (9,10). Fig. shows stress-strain data for 
the first cycle loading, together with the Ramberg-Osgood curve fitted 
these results 0.2, 0.02 percent), and used the present 
study. These Ramberg-Osgood parameters are slightly different from those 
presented Ref. for the same soil. The parameters Ref. represent 
further refinement the interpretation the original test results. However, 
such small variations these parameters have only minor effects the results. 
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The damping provided the Ramberg-Osgood parameters very small 
strains less than about provide damping ratio three percent 
small strains, estimated the basis other laboratory tests soft clays 


reported the literature (19), shear viscosity, was selected for each layer, 
using the equation: 


which: shear viscosity the layer, kips per square feet per second; 
desired damping ratio; 0.03 for all layers; G,,,, shear modulus small 
strains layer, kips per square feet; and period soil profile 
very small strains, estimated equal 0.95 sec for the profile. The values 
the Ramberg-Osgood parameters and the shear viscosities used for all layers 
are summarized Fig. controlled strain cyclic tests the secant modulus 
reduces with number cycles. The ratio the secant modulus Nth cycle 


Nn 


KSF 


DATA 
0.5 


CYCLIC DEVIATORIC STRESSES, (+ Gg),KSF 


os 


CYCLIC AXIAL STRAIN, (+€), PERCENT 
FIG. 3.—Backbone Curve for San Francisco Bay Mud (Ref. 


and first cycle measare degradation the stress-strain relationship for 


soil. The ratio, defined degradation index and related number 
cycles as: 


The degradation parameter, dependent the cyclic strain and has been 
discussed detail Idriss, al. (10). For the sand and stiff very stiff 
clay layers, was assumed the outset that the values did not degrade 
during this earthquake, due the anticipated small strains these layers. 
Degradation the soft clay layer was allowed the analysis, and values 
the degradation parameter for the soft clay were specified based the 
experimental results for soft Bay Mud reported elsewhere (10). 


Site Response 


The response the Southern Pacific Building site was calculated using both 
nonlinear and equivalent linear models. Identical profiles, soil properties and 


GT9 SOFT CLAY SITES 1205 


input rock motions were used both programs DCHARM and SHAKE. The 
input motion used was the S80E component recorded the Golden Gate rock 
site, with the recorded accelerations multiplied factor 0.65 account 
for the different distances the two sites the source (7). Following the 
current state-of-the-art, the motion was input rock outcrop for the SHAKE 
analysis. For the DCHARM analysis, the input rock motion was assumed 
act the base the soil profile. 

Comparisons between recorded and computed acceleration and velocity re- 
sponse spectra near the ground surface using the SHAKE and DCHARM analyses 
are included Figs. and respectively, for spectral damping ratio 
2.5%. For this low damping, the differences between different models are 
accentuated. Examination Figs. and indicates that significantly better 
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FIG. 4.—Computed and Recorded Acceleration Spectra, Southern Pacific Building 


agreement with recorded motions was obtained with the nonlinear compared 
the equivalent linear method. The equivalent linear analysis somewhat 
overestimates the response the low period range, and significantly underesti- 
mates the recorded peak period about sec. However, both the analyses 
predict spurious peak about 0.4 sec. 

For the Southern Pacific Building site, the nonlinear analysis simulates the 
recorded motions better than the equivalent linear method. The improved 
agreement about sec especially significant, because comparison between 
recorded rock (Fig. and soil (Fig. spectra reveals that the peak this 
period represents the main difference attributable the presence the soil. 
The small strain fundamental period the soil profile the Southern Pacific 
Building site 0.95 sec, which increased the SHAKE runs 1.05 sec due 
nonlinear effects; therefore, this peak the spectrum can associated 
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with the vibration the first mode the soil profile. must noted that 
the fundamental period the building also about 1.1 sec and the 
building suffered substantial partition damage, presumably due this peak 
the soil spectrum (23). Due the practical importance this peak, has 
been considered repeatedly the literature (7,23), and the improved prediction 
the nonlinear analysis important result this study. 

The calculated degradation index, the soft layer the end the shaking 
was 0.99, indicating degradation about modulus values the soil 
during this earthquake. Therefore, the improved agreement between computed 
and recorded motions, previously mentioned, due the use nonlinear 
modeling and not the effect degradation. 
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FIG. and Recorded Velocity Spectra, Southern Pacific Building 


For the nonlinear analyses previously considered, the recorded rock motion 
used for the input was measured along the S80E direction, and the computed 
motion was compared with accelerogram recorded soil along the N45W 
direction. These are 35° apart. Therefore, was decided include the other 
recorded horizontal soil component the Southern Pacific Building (N45E) 
the comparisons. Figs. and show comparisons between the computed 
acceleration and velocity response spectra, respectively, and the ranges defined 
the two components the recorded soil motions, for spectral damping 
2.5%. Examination Figs. and suggests that the good agreement provided 
the nonlinear analysis still valid the two recorded components the 
soil motion are considered. 
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where a(t) the acceleration time history. This integral measure the 
energy the accelerogram. Husid (6) plotted the buildup this energy with 
time study the evolution the level shaking. 

addition the comparisons considered earlier between computed and 
recorded response spectra (frequency domain), comparisons the time domain 
were performed. These included comparing the recorded N45W accelerogram 
the Southern Pacific Building site, together with its corresponding Husid 
plot with the corresponding accelerograms and Husid plots calculated the 
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FIG. and Recorded Velocity Spectra, Southern Pacific Building 


nonlinear and equivalent linear analyses. These comparisons are presented 
Fig. The general appearance the calculated and recorded motions are 
very similar, including the existence both strong, high frequency part 
which ends sec-5 sec, followed long period part having smaller 
accelerations and predominant period the oscillations about sec-1.2 
sec. The build-up energy the two records, shown the Husid plots 
the upper part Figs. 8(a) and 8(b), are also very similar. There 
almost straight-line steep increase energy both records until about 
sec, where there break both plots that coincides with the end the 
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strong part the record; after this break, the plots are much flatter, due 
the smaller accelerations the long period part. The total energy the end 
the peak acceleration predicted the nonlinear analysis very similar 
the recorded value (about 0.05 g). The long period part after sec which 
present both the recorded motion and the nonlinear simulation, totally 
absent the equivalent linear simulation [Fig. 8(c)]. The Husid plot from 
the SHAKE analysis Fig. 8(c) completely flat after about sec with 
smaller than that for the recorded motion Fig. 8(a). 

The comparisons Fig. are very useful understanding the origin 
the strong part and the long period part contained the motion recorded 
the site. The accelerations arriving the rock underlying the site became 
very small after about sec-4 sec. This clearly observed the Golden 
Gate rock record shown Fig. Before sec, the response the soil 
strongly affected its nonlinear stress-strain behavior, which caused the soil 
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profile vibrate system with relatively high damping. result, the 
peak input rock acceleration 0.05 was neither amplified nor deamplified 
the soil. After sec, the rock shaking decreased significantly and the soil 
responded almost elastically with low damping. The coincidence between long 
period soil recorded motions and the period the profile (about one second 
both cases) suggests that these rock motions after sec excited the fundamental 
mode the soil profile. This provides reasonable explanation for the peak 
the recorded soil spectrum sec and also for the good agreement between 
Figs. and the other hand, the equivalent linear approach, the 
moduli and damping are chosen compatible with the average effective 
strain. For every cycle strain lower than this effective strain, the damping 
used too high. Since the long period rock motions after sec were 
strain amplitude, they were overdamped the equivalent linear analysis. 
result, the equivalent linear simulation underestimated the soil motions after 


sec; this reflected Fig. 8(c) and the long period ordinates Figs. 4(b) 
and 
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may noted that this apparent deficiency the equivalent-linear procedure 
can overcome also conducting progressive analysis, similar that 
previously used the Union Bay site (20). Strain compatible modulus and 
damping values are applicable for only part the input and are changed 
and updated the input motion progressively applied the analysis profile. 
this way, possible, for example, use higher modulus and lower damping 
values during the latter part the input motion. 

the equivalent-linear analysis the rock motion was input outcrop, 
but the nonlinear analysis the rock motion was input the base the 
soil profile. quantify the importance this effect, additional SHAKE 
analysis was performed, with the input motion the base the soil profile. 
comparison spectral ordinates the recorded accelerations the Southern 
Pacific site with those computed the different analyses shown Table 
for spectral damping 2.5%. When the rock motions are input the 
base the profile the SHAKE analysis instead being input rock 
outcrop, the agreement with the recorded spectra improves the period 


TABLE 1.—Results Analyses for Southern Pacific Building Site 


SHAKE SHAKE 
(1) (2) (4) 


(5) 


Peak spectrum range 

0.2 sec-0.3 sec 0.226 0.223 
Peak spectrum range 

0.9 sec 0.137 0.089 


“SHAKE analysis with rock motion input rock outcrop. 
analysis with rock motion input base soil profile. 


the soil (about 1.1 sec), but deteriorates low periods. Therefore, only 
the nonlinear analysis performed using DCHARM provides, for this case, 
consistently good agreement for wide range periods. 


Station No. 119 Site, 


The city Tokyo, Japan, was shaken the Higashi-Matsugama earthquake 
July 1968. The magnitude the earthquake was about 6.4, its epicenter 
was located about the northwest Tokyo, and had focal depth 
about km. The earthquake was widely recorded the strong motion 
network the Tokyo region (15,16). Information the recorded motions, 
including digitized printouts the accelerograms, response spectra and soil 
information from borings some the stations has been published the 
Association for Science Documents Information, Tokyo (1). The record for 
soft clay selected for this study the component Station No. 119, obtained 
the basement two-story building. The record selected control motion 
the component Station No. 118, obtained the basement five-story 


Computed 
0.063 0.083 
0.065 0.127 
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building dense Tokyo gravel. Both stations were similar distances from 
the epicenter (40 km), and the distance between the sites about 
km. 


Unlike the case the Southern Pacific site studies described earlier, Station 
No. 118 was not rock but shallow stiff soil, and was also some 
distance from the soft site Station No. 119. Therefore, somewhat different 
approach was utilized for the analyses the Station No. 119 site. estimate 
the peak bedrock acceleration the site, available correlations were used for 
magnitude 6.4 event hypocentral distance approx km. value 
0.06 was estimated and was used for the analyses (0.04 had been recorded 
the stiff site Station No. 118). value about 0.35 sec was estimated 
for the predominant period the rock motions. 

Two input rock motions were used the analyses. One the input rock 
motions used was based the component recorded Station No. 118; 
the record was deconvolved using the program SHAKE the Lower Tokyo 
Formation, and the deconvolved accelerogram was then adjusted give peak 
acceleration 0.06 and predominant period 0.35 sec. Acceleration spectra 
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FIG. 9.—Recorded Acceleration Spectra Stiff Site, Tokyo 


the recorded motion Station No. 118 are shown Fig. for spectral 
damping values 2.5% and 10%. The other accelerogram used input rock 
motions was the N21E component the Castaic strong-motion record obtained 
rock during the 1971 San Fernando earthquake, which had magnitude 
about equal that the Higashi-Matsugama event. This record was adjusted 
give peak acceleration 0.06 while the original predominant period 
(0.32 sec) was left unchanged. 

The use two input rock motions, together with the variation soil properties 
examined the following, was designed account for some the uncertainties 


involved the calculations. this way, limited parametric study was performed 
site No. 119. 


The soil profile Tokyo Station No. 119 shown Fig. 10. The profile 
down depth 140 was based boring information given Ref. 
The profile below depth about 140 was based data presented 
Refs. and 21. Bedrock was assumed depth 500 ft. The site overlain 


0.05 ——+——— — 
Spectro! Domping « 100% 
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about silty sand and silty clay, followed 75-ft layer soft 
silty clay. Below the soft clay there medium stiff clay depth 140 
ft, followed 30-ft layer Tokyo gravel. Between 170 and the tertiary 
bedrock there geologic unit known Lower Tokyo Formation. 

The soil and rock properties used the analyses Station No. 119 are 
listed Fig. 10. The values unit weight, y,, and maximum shear modulus, 
estimated based literature survey Tokyo subsoil data and 
correlations established Ohsaki (14). The Ramberg-Osgood parameters 
were arrived similar way that described for the Southern Pacific 
Building site profile. For granular and stiff very stiff materials, Ramberg-Osgood 
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FIG. 10.—Soil Profile and Parameters for Response Analyses, Tokyo Site No. 119 


parameters 1.0, and 0.02% were selected. For the medium 
stiff clay, and 0.02% were used. The parameters 
used for the soft silty clay were the same for San Francisco Bay Mud 
the Southern Pacific Building site analyses. Only the soft silty clay was 
allowed degrade. The modulus degradation parameter for the soft layer 
was assumed the same that used for the San Francisco Bay Mud. 
shear viscosity was selected for all layers provide damping ratio 
3%. Due the uncertainties involved the selection soil parameters for 
the analyses, two sets values were used define 
and Bound”’ stiffnesses, and these are shown Fig. 10. For SHAKE 


a Medium Stiff Silty Cloy 
140 
170 
17,000 
i 
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TABLE 2.—Station No. 119 Site, Tokyo Cases Studied 


Input rock motion 


Soil 
(2) (3) 


Best estimate 
Upper bound 
Best estimate 
Upper bound 


Site no. 118 (deconvolved) 
Site no. 118 (deconvolved) 
Castaic 
Castaic 


Fig. 10. 
TABLE 3.—Comparison Recorded and Calculated Peak Accelerations 


Nonlinear 
response 
calculations, 
gravitational gravitational gravitational 
acceleration acceleration acceleration 


(2) (3) (4) 


Equivalent 
linear 


Recorded, method, 


Southern Pacific site 
Tokyo site no. 119 
case 


Recorded (NS) 


Nonlineor Anolysis 
Case B (See Tobie |) 


z 
z 
oO 
=) 
fo} 
a 
< 


1.00 
Recorded (NS) 


Equivalent Linear Analysis 
Cose 8 (See Tobie |) 


ABSOLUTE ACCELERATION IN G's 


0.10 0.30 1.00 3.00 
PERIOD IN SECONDS 


10.00 


FIG. and Recorded Acceleration Spectra, Tokyo Site No. 119 


Case 
(1) 
Site 
(1) 
0.047 0.054 0.062 
0.050 0.035 0.050 
0.050 0.038 0.064 
0.050 0.024 0.036 
0.050 0.034 0.050 
0.25 
af \ Spectral Damping 
0.03 3.00 10.00 
0.30 
A \ 
\ 
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runs the best estimate with the modulus and damping curves 
Seed and Idriss (19) were used. 


Response 


The response the Tokyo Station No. 119 site was calculated using the 
nonlinear and equivalent linear methods analysis (programs DCHARM and 
SHAKE, respectively). Identical profiles and soil properties were used both 
sets calculations. Two input rock motions were assumed act the base 
the soil profile, and two soil stiffnesses were used for the analyses, previously 
described. Therefore, the four cases A-D, described Table were analyzed 
using both the nonlinear and the equivalent linear methods. 

Comparison between the computed and recorded response spectra near ground 
surface are shown Figs. 11-16. the Southern Pacific case, spectral 
damping value 2.5% was used for all comparisons. Figs. and show 
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comparisons between the recorded and computed acceleration and velocity 
response spectra for Case Examination these spectra indicates again that 
the nonlinear analysis provides improved agreement with the recorded values. 
This true especially for periods less than 0.8 sec, where the equivalent linear 
method overestimates the response. 

Figs. and present the acceleration response spectra for Cases A-D 
using the equivalent linear and nonlinear analyses, respectively. Fig. compares 
the ranges results corresponding the two analytical models obtained from 
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FIG. 14.—Computed Acceleration Spectra, Nonlinear Analyses, Tokyo Site No. 119 


Figs. and 14. Fig. compares the acceleration spectrum recorded Station 
No. 119 with the range computed nonlinear analysis. The recorded motion 
plots inside the computed range very close all periods longer than 
0.2 sec. examination Fig. indicates that the recorded spectrum coincides 
approximately with the upper bound the analytical range for periods longer 
than about 0.8 sec. This upper bound defined the two analyses using 
rock input the deconvolved Station No. 118 motion, while the computed 
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spectra using Castaic rock input are significantly lower that period range 
(see Fig. 14). 

These results indicate that response analyses should preferably used within 
the framework parametric study. only one analysis performed 
one rock motion assumed, the calculated results may unrealistic. 

comparison between recorded spectra the Station No. 118 stiff site (Fig. 
and the soft Station No. 119 site (Fig. 16) suggests that the main differences 
response attributable the presence the soil are the period range 
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FIG. 15.—Comparison Acceleration Spectra Computed Nonlinear and Equivalent 
Linear Analyses, Tokyo Site No. 119 
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FIG. 16.—Acceleration Response Spectra Recorded and Computed Nonlinear 
Analyses, Tokyo Site No. 119 


between approx sec and sec. The fundamental period the soft site, obtained 
the four SHAKE runs and including nonlinear effects, was computed 
the range 1.2 sec-1.8 sec, with higher values for Cases and 
(Best Estimate stiffness) than for Cases and (Upper Bound stiffness). 
Therefore, the Southern Pacific example, the ground motions the soft 
site the sec period range seem have been amplified the 
characteristics the soil profile. interesting, however, that the computed 
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response this period range seems more sensitive changes the 
input motion than variations the soil properties. For all four nonlinear 
analyses, the calculated value the degradation index the soft layer 
the end shaking was 0.99, indicating negligible degradation modulus 
values the soil during this earthquake. 


Review ano 


The main purpose this study was calibrate recently developed nonlinear 
soil response procedure (9,10). the two sites selected for this purpose, 
strong-motion records having peak accelerations about 0.05 were available 
soil profiles containing layer soft clay, and simultaneous records 
rock stiff sites were also available, which are used control motions. One 
important difference between the two cases was that for the Southern Pacific 
Building site, the input motion and subsurface soil conditions and properties 
were well known, whereas, the Tokyo sites, they were estimated using 
correlations. 

The comparisons with the recorded motions the two soft sites suggest 
that the use DCHARM provides improved estimates the earthquake ground 
response, and therefore improvement over the current use SHAKE. 
The results the analyses indicate that the equivalent linear method computes 
systematically higher spectral ordinates than the nonlinear method low periods. 
This may due the fact that, program DCHARM, 0.03 for all 
layers. The shear viscosity may responsible for overdamping the high frequency 
response. The analyses conducted part this study show reasonable 
value assume for overall response results. For the long period range, the 
nonlinear method gives response which higher about equal that computed 
the equivalent linear calculations. the case the Southern Pacific site, 
the response obtained the long period range using the nonlinear method was 
reasonable compared the recorded values and was related the response 
the soil after about sec shaking. 

good agreement was obtained between ground response computed the 
nonlinear method, and the recorded ground motions the Southern Pacific 
Building site. This agreement was good both the frequency domain (response 
spectra) and time domain, including the simulation the strong part followed 
the long period part, which appear the recorded accelerogram. This good 
agreement may related the availability both realistic rock motion input 
and measured soil parameters. the Station No. 119 site Tokyo, reasonably 
good agreement was also obtained between computed and recorded ground 
motions some the nonlinear analyses. However, parametric study, including 
possible variations assumed rock motions and soil properties, showed that 
the quality the agreement was not good for the Southern Pacific Building 
site. This emphasizes the need for parametric studies when using response 
calculations engineering applications. 

all nonlinear analyses performed this study, there was practically 
modulus degradation the soft clay layer for the level acceleration and 
duration earthquakes studied. This suggests that substantially higher accelera- 
tions and durations are needed cause significant degradation soft clay 
deposits. Additional studies are needed evaluate the levels earthquake 
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shaking intensities which degradation becomes important factor the 
behavior soft clay sites. believed that nonlinear analyses such cases 
will yield results closer reality than the equivalent linear methods. 
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ANISOTROPIC STRENGTH COHESIONLESS SANDS 
Masanobu 


INTRODUCTION 


Natural situ soil deposits are neither isotropic nor homogeneous. They 
often possess layered structures, e.g., alternating layers having quite different 
particle size distribution, where each homogeneous layer nonspherical particles 
have strong tendency exist parallel nearly parallel the layered structures 
(3,9,20). The anisotropy and heterogeneity may have considerable influence 
the strength and other properties soils; e.g., the shear strength clays 
(6,22,23), and sands (2,3,11,15,20,21), and the bearing capacity sand 
(12,20,22). 

this paper, the writer reports some results plane strain tests samples 
having layered structures well being anisotropic. The results may provide 
guidance for estimating the internal friction angle soils with complicated 
properties. 


Following Brewer (5), fabric mean the spatial arrangement solid 
particles and the associated voids. Oda (15) suggests that the fabric sand 
can defined the following two elements: (1) Orientation nonspherical 
particles (orientation fabric); and (2) mutual relation particles (packing). 

Orientation Fabric.—This means the parallel alinement nonspherical partic- 
les. The thin section method (11,13,15) and X-ray method (1,2) were proposed 
order obtain quantitative estimates the orientation fabric. 

Packing.—Suppose that particle has contact with four particles 
four contacts reasonable expect that the relation with 
the neighboring particles G,-G, given by: (1) The four unit vectors, n,-n,, 
which are normal the contact areas C,-C,; and (2) the number contacts 
per particle, Since very large number particles involved, introduce 
density functions, E(n) and f(m), define the statistical distribution and 
the considered soil mass. The probability density function, showing 
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the distribution the unit normals has been considered Field (7), Oda 
(15,16,17), Oda and Konishi (18), and Mulilis, al. (13). The probability 
distribution the number contacts per particle, has been used 
Smith, al. (25), Bernal and Mason (4), Field (7), Oda (17), O’Conner, 
al. (14), and Grivas (8). 

From the detailed observation the fabric sand, Oda (15) showed that 
two types anisotropy can distinguished uniform sand which deposited 
under the action gravitational force; these are: (1) Anisotropy the preferential 
alinement nonspherical particles parallel the horizontal; and (2) anisotropy 
the concentration the unit vectors, the vertical direction. Both types 
anisotropy are possible naturally deposited sand nonspherical particles, 
while the mass spheres only anisotropy due the concentration n’s 
only possible. 

include heterogeneous mass our consideration, one more type 
anisotropy must added; i.e., (3) anisotropy the layer structure (3) 
which formed the alternative deposition different materials such 
sand, silt, and clay. Since the three types anisotropy are usually interrelated 
complicated manner, great importance consider separately the 
effect each anisotropy the overall mechanical anisotropy soils. 


Test AND SPECIMENS 


Test Materials.—Test materials are Toyoura sand, Tochigi sand, and glass 
beads, the physical properties which are summarized Table [An axial 


TABLE Properties Test Materials 


Axial 
Void Ratios Ratio 


Grain size, 
millimeters 


(6) 


Test materials 


(1) 


Roundness 


(5) 


Toyoura sand 0.99 0.63 0.61 0.42 0.105 
Tochigi sand 0.42 0.25 
Glass beads 0.42 0.105 


ratio, Table introduced identify the average shape particles 
(15)]. sand which characterized high value, close unity, consists 
particles close spheres. According the previous study natural sands 
(21), the corresponding R-values range from 0.5-0.7. 

Toyoura sand uniform sand with high content feldspar mineral. 
This consists rod-like particles with subangular surfaces. Tochigi sand 
also uniform sand, but consists quartz grains with subangular surfaces. 
This selected testing material order compare its results with those 
Toyoura sand, because both have similar physical properties, but different 
anisotropy. The sample glass beads consists uniform glass balls. This 


artificial material used study the contribution n-concentration the 
strength anisotropy. 
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Description Samples.—Samples used are classified the basis their 
uniformity and anisotropy, follows. 

Uniform-Anisotropic Samples.—Oven dried sands are poured into sample 
former through funnel which set 900 above the sample former. 
This pouring method yields uniform-dense sample with relative density 
90%. clear from the published data (15,21) that the resulting fabric 
Toyoura sand, with the void ratio highly anisotropic not only 
the preferred alinement rod-like particles parallel the horizontal but 
also the high concentration the vertical. The sample Tochigi 
sand, with the void ratio does not show appreciable parallel 


INITIAL 
HORIZONTAL 
PLANE 


GLASS 
BEADS 


SAND 


. 


FIG. 1.—Ratio, between Thickness Toyoura Sand Layer and that Glass Beads 
Layered Samples, and Inclination Angle, Initial Horizontal Plane Maximum 
Principal Axis 


alinement nonspherical particles, but shows the same concentration 
the vertical that Toyoura sand. Kallstenius and Bergau (10) first reported 
the anisotropic fabric random mass complete spheres, stating that 
mass (of complete spheres) anisotropic the extent that the height each 
tetrahedron less than required for This finding corresponds 
the fabric anisotropy the concentration n’s mass spheres. 

summary, samples made the pouring procedure are characterized 
the uniform-anisotropic fabric irrespective the shape particles. 

Uniform-pseudoisotropic samples: Oven-dried Toyoura sand poured make 
sand layer thickness. After each pouring, hand plunger with 
diam plunged directly into the sand destroy the initial anisotropic 


5 
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arrangement particles. Each penetration carefully controlled less 
than depth, and total 100 penetrations—200 penetrations made 
per each layer. The side wall the sample former tapped order obtain 
the same void ratio, ranging from that the sample made 
the pouring method. This procedure repeated until the sample formation 
completed. can say, the basis the fabric analysis Oda (15), that 
very difficult make ideally isotropic samples destroying the pre-existing 


kgf/cm? 


FIG. 2.—Stress-Strain Relations Uniform-Anisotropic Toyoura Sand 


anisotropic fabric. The samples, however, are much closer being isotropic, 
especially with respect the distribution n’s, than the samples which are 
prepared the pouring method. 

Layered-Anisotropic Samples.—Toyoura sand and the glass beads are deposited 
alternately using the pouring method. Accordingly, the samples this series 
consist successive layers uniform-anisotropic Toyoura sand and glass beads 
(Fig. 1). The thickness ratio, defined 
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thickness the glass beads 


thickness Toyoura sand 


controlled 0.1 1.0 changing the thickness the glass-bead layers. 
The thickness Toyoura sand kept constant mm. 

Layered-Pseudoisotropic Samples.—These samples consist successive layers 
uniform-pseudoisotropic Toyoura sand and glass beads. Each successive layer 
pseudoisotropic because its initial anisotropy destroyed the plunging-tap- 
ping method. All samples had 1.0 for this series. 


TOYOURA SAND 
TOCHIGI SAND 


2 
0.5 


FIG. 3.—Anisotropic Strength Response Homogeneous-Anisotropic Samples 


Plane Strain Test.—Drained plane strain tests are performed 
40.4 samples previously mentioned. The direction the maximum 
principal stress, inclined the principal plane anisotropy (parallel 
the initial horizontal plane) angle, order determine the effect 
the anisotropic fabric and the layer structure the anisotropic shear strength 
(Fig. 1). The detailed description for making the samples with different inclination 
angles, and for experimental procedure has been given Oda, Koishikawa 
and Higuchi (21). 

The inclination angle, selected 0°, 15°, 22°, 30°, 45°, 60°, 75°, 
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and 90°, and the minimum principal stress, kept constant during each test, 


Uniform-Anisotropic Samples.—Fig. shows typical stress-strain relations 
the uniform-anisotropic Toyoura sand (void ratio being 0.659-0.663). 

rather difficult choose unique coefficient which the shear strength 
anisotropic sand can fully defined. friction angle, defined 


used shear strength index, which the suffix refers the stress 
failure. 

The experimental results for the uniform-anisotropic samples Toyoura sand, 
Tochigi sand and glass beads are summarized Fig. plotting the friction 
angle, against the inclination angle, the principal plane the fabric 
anisotropy (parallel the initial horizontal plane) the maximum principal 
stress axis. The samples belonging the same series have almost the same 
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FIG. 4.—Stress-Strain Relations Uniform-Pseudoisotropic Toyoura Sand 
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void ratio; Toyoura sand 0.66 0.67, Tochigi sand 0.74 0.77, glass 
beads 0.58 0.61. can say, therefore, that the corresponding change 
with respect the angle shows the effect the anisotropic fabric. 

The most important observation Fig. the remarkable anisotropic response 
the series the glass beads, consisting almost complete spheres. The 
change due the anisotropy much that Toyoura sand which 
consists rod-like particles with their strong alinement parallel the initial 
horizontal plane. Based the analysis experimental results from conventional 
triaxial tests, Oda (15) concluded that the anisotropic distribution nis crucial 


~— 


PLUNGING PER LAYER 


ANISOTROPIC 
TOYOURA SAND 


FIG. 5.—Pseudoisotropic Strength Response Homogeneous-Pseudoisotropic Sam- 
ples 


the determination the anisotropic response. This idea seems supported 
the fact that glass beads with orientation fabric show anisotropic 
shear strength similar that Toyoura sand with strong orientation fabric. 
the same time, can also point out that Toyoura sand with 0.61 
more anisotropic its strength behavior than Tochigi sand with 0.71. 
explain the difference from microscopic point view, cannot neglect 
the effect the orientation fabric the anisotropic response. Although the 
orientation fabric only one the influential factors controlling the function 


30° 60° 90° 
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nonspherical particles than mass spherical particles. 

Another important observation Fig. that each series has minimum 
angle ranging from 15°-30°, depending the material property. 
should noted that such minimum was observed with conventional triaxial 
compression tests (2,3,15). Particles can move easily without any strain confine- 
ment the triaxial compression tests, while movement particles restricted 
the plane strain tests, because zero strain must maintained the direction 
the intermediate principal stress, This also responsible for the fact 
that the anisotropic response more pronounced the plane strain condition. 
Because the easy mobility particles the triaxial compression test, the 


60° 
30° 


FIG. 6.—Stress-Strain Relations Layered-Anisotropic Sample 


initial anisotropy lost during the course deformation and final fabric 
failure attained irrespective the initial fabric. the other hand, the 
initial anisotropic fabric well preserved failure the plane strain tests. 
This means that the principal plane the initial anisotropy weakness plane 
that the sample with 15° 30° more easily fails along this plane. 
Uniform-Pseudoisotropic Samples.—Fig. shows typical stress-strain relations 
the pseudoisotropic Toyoura sand, which plunged 200 times per layer 
destroy anisotropic fabric. Samples quoted Fig. have void ratios ranging 
from which almost correspond the void ratio the uniform-ani- 
sotropic Toyoura sand (Fig. 2). clear, Fig. compared with Fig. 
that the strength anisotropy completely destroyed the plunging procedure. 
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Results this series are summarized Fig. The open circles are obtained 
from the samples which are plunged 200 times per layer, while the closed circles 
are obtained from the samples which are plunged 100 times per layer. broken 
line Fig. shows the experimental result the uniform-anisotropic Toyoura 
sand. 

Since the angle depends the angle with minimum value 
60°, cannot say that the sample ideally isotropic. The plunging-tapping 
method does not yield ideally isotropic sample. However, clear that 
the dependence much less the uniform-pseudoisotropic sample 
than the uniform-anisotropic sample. Note that the mean value 


~ 
° 


FIG. 7.—Anisotropic Strength Response Layered-Anisotropic Samples 


the uniform-anisotropic Toyoura sand almost equal the angle the 
Toyoura sand with 90°. This result suggestive 
for the determination the mean value applicable anisotropic 
sand. 

Layered-Anisotropic Samples.—Fig. shows typical stress-strain relations 
the layered-anisotropic samples with 1.0. Fig. also shows the results 
for the layered-anisotropic samples with 0.1 and 1.0, plotting 
against These data Fig. are located within the zone bounded two 
curves. The upper broken curve corresponds the result the uniform-aniso- 
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tropic Toyoura sand, while the lower broken curve corresponds the result 
the uniform-anisotropic glass beads. 

The strength the layered-anisotropic samples decreases with the increase 
the content glass beads. The decrease, however, seems not depend 
linearly the ratio (3). Note that the drastic decrease clearly observed 
the sample with 15° 30°. 

Layered-Pseudoisotropic Samples.—Tests are performed only samples with 
1.0, shown Fig. dotted line the same figure the result 
the layered-anisotropic samples. The observed anisotropy the layered-pseu- 
doisotropic samples corresponds the layer structure, while that the layered- 


LAYERED- PSEUDOISOTROPIC 
SAMPLE 


FIG. 8.—Anisotropic Strength Response Layered-Pseudoisotropic Samples 


anisotropic samples corresponds the combined effect the fabric anisotropy 
each layer and the layer structure. Both results, however, are similar, except 
each corresponding minimum different value The similarity 
the results seems support the idea that the layer structure predominant 
factor for controlling the strength anisotropy. 


Natural sand deposits are generally anisotropic well heterogeneous. 
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order provide insight regarding fundamental factors affecting the strength 
behavior such complex soils, have performed plane strain tests 
uniform-anisotropic, uniform-pseudoisotropic, layered-anisotropic and layered- 


pseudoisotropic samples. The experimental results support the following conclu- 
sions: 


can distinguish three forms fabric anisotropy; e.g., orientation, 
packing, and layering. Fabric anisotropy due packing crucial the 
determination strength anisotropy homogeneous sand. Fabric anisotropy 
due parallel alinement particles only one the influential factors 
controlling the packing. Layering can another predominant factor the 
determination strength anisotropy. 

should expect anisotropic strength response even uniform mass 
composed spheres. This caused the anisotropic distribution contact 
normals which can take place even random mass spheres. 

necessary always think that assembly particles generally anisotropic. 

Tests pseudoisotropic samples whose anisotropic fabric destroyed 
applying plunger give reasonable estimates the mean strength for 
anisotropic sands. want calculate the bearing capacity uniform-aniso- 
tropic foundations the conventional way (26), e.g., testing the pseudoiso- 
tropic sample may advisable for obtaining reasonable estimate the internal 
friction angle. 
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The following symbols are used this paper: 


void ratio; 
E(n) probability density contact normal 
unit vector normal contact; 
axial ratio showing particle shape; 
= 


inclination angle initial horizontal plane; 
maximum principal strain; 
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volumetric strain; 

thickness ratio; 

maximum principal stress (in effective); 
intermediate principal stress (in effective); 
minimum principal stress (in effective); and 
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GEOTEXTILE-REINFORCED UNPAVED 
Jean-Pierre ASCE, and Laure 


Geotextiles the name now universally adopted for fabrics used geotechnical 
engineering. They have been used extensively unpaved roads. Geotextiles 
are placed between the aggregate layer and the subgrade soil (Fig. perform 
several functions: (1) Separation; (2) filtration; (3) drainage; and (4) reinforcement. 

providing reinforcement, geotextiles improve the performance unpaved 
roads: for given thickness aggregate layer, the traffic can increased 
or, other words, for the same traffic, the thickness aggregate layer can 
reduced (in comparison with the required thickness when geotextile 
used). 

This improvement performance has been recognized users the field, 
but design procedures are lacking. The study presented this paper addresses 
this lack developing method which enables the engineer calculate the 
required thickness aggregate layer and make the proper selection the 
geotextile used. The results are presented the form charts established 
using combination of: (1) Formulas relating aggregate thickness and traffic 
for unpaved roads without geotextile; and (2) quasi-static analysis comparing 
unpaved roads behavior with and without geotextile. 

The readers and possible users these charts must recognize that the present 
study applies only purely cohesive subgrade soils and mostly applicable 
roads subjected light medium traffic truck passages over 
the lifetime the road). 


Derinitions 


Geometry Unpaved Road.—The cross section unpaved road defined 
Fig. which: thickness aggregate layer when there geotextile 
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(in meters); thickness aggregate layer when there geotextile (in 
meters); and reduction aggregate thickness resulting from 
the use geotextile (in meters). 

The subgrade soil assumed homogeneous, least over thickness, 
Sufficient for the development plastic zone (see later, Eq. 19). 

Definition Axle Load and Contact Area.—Any vehicle load can replaced 
equivalent single axle load using the procedure described Appendix 
Therefore, the analysis presented this paper refers single axle. The 
geometry this axle depicted Fig. 2(a) which distance between 
the midpoints the two sets wheels (in meters). Dual wheels are considered 
because they are more common than single wheels for trucks using unpaved 
roads. 


The axle load, considered evenly distributed between the four 
wheels: 


which: axle load (N); contact area tire (in square meters); 
and tire inflation pressure (in pascals); assumed equal the average 


FIG. 1.—Geometry Unpaved Road FIG. 2.—Geometry Vehicle Axle and 
Definition Contact Area: (a) Geo- 
metry Vehicle Axle; (b) Tire Contact 
Areas; and (c) Equivalent Contact Area 
used Analysis 


value the actual contact pressure (nonuniformly distributed) between each 
tire and the aggregate layer. 

The soil between the tires dual wheel mechanically associated with 
these tires [Fig. 2(b)]. Since failure the aggregate layer and subgrade 
soil can occur between the tires, each double contact area 2A. replaced 
the theoretical study rectangle larger area [Fig. 2(c)]. 
examining several dual tire prints, the following value appears reasonable: 


The actual contact pressure (nonuniformly distributed) between each tire and 
the aggregate layer induces the same mechanical effects the subgrade 


contact (assumed uniformly distributed) between 
rectangle and aggregate layer; therefore: 


GT9 


subgrade soi! 
(2) 
(3) 
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The relationship between equivalent contact pressure, p.., and tire inflation 
pressure, deduced from Eq. and 


Because the subsequent analysis bidimensional, must eliminated. 
Examination typical dual tire prints leads the following approximate value 
for on-highway trucks: 


the case off-highway trucks, the tires are very wide and the following 
value more appropriate: 


for on-highway trucks 


for off-highway trucks 


c 


Eqs. and are useful for the subsequent analysis because vehicles are 
usually characterized axle load and tire inflation pressure, p., which 
almost equal the actual contact pressure previously stated. 

Properties Aggregate.—The aggregate assumed have the properties 
usually required ensure proper distribution the applied load. Practically, 
this means that the California Bearing Ratio (CBR) the aggregate layer 
larger than 80. 

Properties Subgrade Soil.—The subgrade soil assumed saturated 
and have low permeability (silt, clay). Therefore, under quick loading (such 
traffic loading), behaves undrained manner. Practically, this means 
that: (1) The subgrade soil incompressible; and its friction angle nil. 
Consequently its shear strength equal its undrained cohesion, c,. 

The value measured the laboratory using unconsolidated undrained 
triaxial test quick shear test, and the field using vane shear test. The 
value can also evaluated using cone penetrometer deduced from 
CBR value using one the following relationships: 


(in Pascals) 30,000 CBR 


which: cone resistance (in Pascals); and CBR California Bearing 
Ratio (dimensionless). These two equations are widely accepted and their 
justification beyond the scope this paper. 


Pec V2 (4) 
Eliminating from Eq. and using Eq. yield: 
(9) 
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Properties geotextile mechanical properties relevant 
this study are best described the geotextile tension-elongation behavior 
biaxial tensile test, where lateral deformation the geotextile sample 
prevented. Tension (in Newtons per meter) defined tensile force (in Newtons) 
per unit width (in meters). 

Fig. presents typical tension-elongation curve and shows the definition 
the secant modulus, the following, only the secant modulus the 
transverse direction the road considered. Typical values for elongations 
between and 20% are: kN/m for plain rubber plastic sheets; 
kN/m for nonwoven geotextiles; and kN/m for 
woven geotextiles. 

The geotextile elongation associated with road deflection must smaller 
than the elongation failure 


(N/m) 


FIG. 3.—Typical Tension-Elongation Curve Geotextile Elongation Failure; 


Point 


For woven geotextile, the tension-elongation curves obtained from uniaxial 
test (as usually performed the textile industry) and biaxial test are 
not very different. For nonwoven geotextile, the difference usually important 
and the biaxial test strongly recommended. 

addition tensile behavior, necessary address the friction charac- 
teristics the geotextile. beyond the scope this paper elaborate 
this subject. This study assumes that the geotextile rough enough prevent 
failure the aggregate layer sliding along the geotextile surface. 


Mechanisms Failure Unpaved Road.—The failure unpaved road 
can occur three different locations: (1) Aggregate layer; (2) foundation soil; 
and (3) geotextile (if any). The present study does not consider failure 
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the aggregate layer. other words, assumed that: (1) The friction coefficient 
the aggregate layer large enough ensure the mechanical stability 
the layer; and (2) the friction angle the geotextile contact with the aggregate 
layer under the wheels large enough prevent sliding the aggregate layer 
the geotextile. 

Basically, this study evaluates the risk failure the foundation soil and 
the geotextile. 

Load Distribution Aggregate Layer.—The aggregate layer assumed 


provide pyramidal distribution with depth the pressure, p.., applied 
its surface (Fig. 4). 


Therefore: 


which: weight per unit volume the aggregate (in Newtons per cubic 
meter); pressure the base aggregate layer when geotextile 


subgrade soil 


FIG. 4.—Load Distribution Aggregate Layer: (a) Case Without Geotextile; and 
(b) Case with Geotextile 


used (in Pascals); pressure the base aggregate layer when geotextile 
used (in Pascals); and and angles load distribution pyramid (in 
degrees) (see Fig. 4). 

From 11, and 12: 


There some uncertainty about the value adopted for and However, 
different approaches lead similar values: (1) experimental approach consists 
considering the shape the deformed aggregate layer after tests (2,4); this 
suggests approximate values tan and tan between 0.5 and 0.7; and 
classical value for the inclination the shear planes the aggregate layer 
were punched; for between 25° and 45°, this yields values tan and 
tan between 0.47 and 0.64. 


aggregate Pec 
geotextile 
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Preliminary calculations have shown that the value the calculated aggregate 
thickness not significantly influenced and long tan and 
tan are between 0.5 and 0.7. Therefore, tan tan 0.6 has been 
assumed all subsequent calculations. 

Bidimensional Approach.—The analysis will bidimensional vertical 
cross section the road. This based upon the following two points: (1) 
The wheels are assumed travel always along the road the same track, 
that every cross section the road receives the same amount loading 
and deforms the same way; thus, there residual deformation the 
direction the road and every cross section has the same mechanical behavior; 
and (2) classical practice soil mechanics evaluate the limit equilibrium 
soil using bidimensional analysis; first, because the tridimensional limit 
equilibrium soil still not well known; secondly, because according 
studies this subject, the difference between bearing capacities deduced from 
bidimensional and tridimensional equilibria frictionless soil not important. 

Elastic and Plastic Behaviors Soil under Loading.—Consider frictionless 
soil subjected uniformly distributed pressure, over limited portion 
its surface and uniformly distributed lateral surcharge, [Fig. 


(M+ de 


efter Whitman 
and Woes (5) 


frictionless soi! 


FIG. 5.—Elastic and Plastic Behaviors Soil under Loading: (a) Loading Pattern; 
and (b) Settlement, Center Versus Applied Pressure 


Classical bidimensional calculations based theory elasticity lead the 


value which induces stresses the soil corresponding the elastic 
limit the soil: 


which: elastic bearing capacity the soil (in Pascals); and undrained 
cohesion the frictionless soil (in Pascals). 


Similar calculations using theory plasticity lead q,,,, maximum value 


which: ultimate bearing capacity the soil (in Pascals). 

Fig. shows the settlement, occurring the center the loaded portion 
the soil surface function the difference, This settlement 
small when the value smaller than and large when the value 

Fig. 6(a) shows the development plastic zones frictionless incompressible 


soil under load equal the ultimate bearing capacity the soil and applied 
over width 2(a). 


UNPAVED ROAD DESIGN 
The depth the plastic zones given by: 


Theory plasticity will used later describe the behavior the subgrade 
soil when geotextile used. this case, according Fig. 


According Eqs. and 18, the minimum value for the subgrade soil thickness, 
ensure the validity subsequent calculations based upon theory plasticity, 
is: 


The width plastic zones, either sides the loaded area, equal 
the width the loaded area. Therefore, two loads are distance smaller 
than [Fig. 6(b)], the heave the soil between the two loads restricted 
and more soil moves outwards than inwards. This consideration important 


FIG. 6.—Mechanism Plastic Failure FIG. 7.—Increase Lateral Surcharge 
Soil: (a) One Load; and Heave Soil 
(b) Two Loads 


will shown later. The width plastic zones also governs the required 
lateral extent the geotextile under the road shoulders. 

conclusion, the soil settles under the load (Fig. and heaves laterally 
(Fig. 6). This phenomenon whole called deflection the soil surface. 
more rigourous analysis would consider that heave the soil increases the 
lateral surcharge (Fig. 7). fact, has been checked that such increase 
does not have significant effect the results the subsequent analyses. 
Therefore, will neglected. 


Quasi-Static aforementioned results concerning theories 
elasticity and plasticity can used the case loaded aggregate layer. 
When there geotextile [Fig. (4a)], the load applied the subgrade soil 
dual wheel and the weight aggregate layer equivalent uniform 
pressure, p,, over width 2h, tan a,), and uniform lateral surcharge, 
elsewhere (with weight per unit volume aggregate, Newtons 
per cubic meters). 
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exceeds the elastic bearing capacity the subgrade soil, the deflection 
will large and complete failure the road will occur after few passages 
the dual wheel. Consequently, the design unpaved road without geotextile 


must consider maximum pressure the subgrade soil equal the elastic 
bearing capacity. Therefore, according Eq. 15: 


will seen later that, the case geotextile-reinforced, unpaved 
road, because different kinematics, the maximum pressure the subgrade 


soil can practically equal the ultimate bearing capacity. 
From and 20: 


(21) 
2(B 2h, 2h, tan 


Therefore, for on-highway trucks, according Eqs. and 


(22) 
Similarly, for off-highway trucks, according Eqs. and 


which: cohesion subgrade soil (in Pascals); axle load (in Newtons); 
tire inflation pressure (in Pascals); thickness aggregate layer (in 
meters); and angle load distribution pyramide (in degrees) (Note: tan 
0.6, previously stated). 

The chart presented Fig. has been established using Eqs. and 23. 
This chart allows determination aggregate thickness, when subgrade soil 
cohesion, and axle load are known. Design engineers should bear mind 
that the use this chart restricted very light traffic (say 1-20 passages) 
because has been established using quasi-static analysis. procedure 
take heavier traffic into account presented the next section. 

Influence Traffic.—An extensive test program unpaved roads (without 
geotextile) has been conducted the Corps Engineers (1). Failure criterion 
was rut depth 0.075 From the test results, Webster and Alford established 
chart (3) giving the thickness aggregate layer function number 
passages and CBR subgrade soil, for standard axle load, P,, 


The following formula has been found good agreement with Webster 
and Alford’s chart: 


0.19 log (24) 


which: thickness aggregate layer (in meters); number passages 


GT9 UNPAVED ROAD DESIGN 1241 


standard axle (load kN); and CBR California Bearing Ratio 
subgrade soil. The use Eq. restricted standard axle load and 0.075 
rut depth. Eq. appears better agreement with Webster and Alford’s 
chart than other proposed formulas giving function log (CBR). 

extend other axle loads, the following equation proposed: 


P 3.95 
i 


which: number passages standard axle (load kN); and 
number passages axle load Eq. (or similar equations) 
widely accepted for traffic paved roads, and because lack any similar 
equation specific unpaved roads, used here. 


Legend 


—-— tire inflation pressure p. = 480 kPa 


—— tire inflation pressure p, = 620 kPa 


= 870 kw 
= 560 kN Eq. 23 
= 390 kw 


fq. 22 


= 80 
| 


(kPa) 


FIG. 8.—Aggregate Thickness Versus Subgrade Soil Cohesion (Quasi-Static Analy- 
sis for Case Without Geotextile) 


log [log 2.34 0.075)]. This expression has been empirically 
deduced from test results presented Webster and Watkins (4), showing that 
increase rut depth with number passages much more marked soon 
rut depth exceeds 0.075 

Combining Eqs. 10, 24, and replacing log the preceding given: 
expression 


0.63 
Cy 
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has been established extrapolation and should not used when the number 
passages exceeds 10,000. 

Fig. presents values obtained from Eq. the case standard axle 
load and 0.3 rut depth. Similar charts have been established for several 
values axle load and rut depth, 

Fig. shows that the curves versus obtained from quasi-static analysis 
(Eq. 22) correspond number passages the order magnitude 
100. With allowable rut depth 0.075 the curves versus would 
correspond number passages the order magnitude 10. This 


nh for tire inflation pressure %° 620 kPa 


hy for tire inflation pressure p. = 480 kPe 


5 


CBR 


FIG. 9.—Aggregate Thickness Versus Subgrade Soil Cohesion (Case without 


Geotextile when Traffic Taken into Account) (Curves Versus are Extracted 
from Fig. 


comparison shows that quasi-static analysis well suited light traffic, 
which seems reasonable. 


Kinematics.—As previously mentioned, the subgrade soil incompressible 
and settlement under the wheels causes heave between and beyond the wheels 
(Fig. 10). Therefore, the geotextile exhibits wavy shape; consequently, 
stretched. 

When stretched flexible material has curved shape, pressure against its 
concave face higher than pressure against its convex face. This known 
Therefore: (1) Between the wheels (BB Fig. 10) 
and also beyond the wheels (AC Fig. 10), although lesser extent, the 
pressure applied the geotextile the subgrade soil higher than the pressure 
applied the aggregate layer the geotextile; and (2) under the wheels (AB 
Fig. 10), the pressure applied the geotextile the subgrade soil smaller 
than the pressure applied the wheels plus the aggregate layer the geotextile. 

the geotextile provides two beneficial effects: (1) Confinement 
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the subgrade soil between and beyond the wheels; and (2) reduction the 
pressure applied the wheels the subgrade soil. These two effects will 
evaluated the following sections. 

Subgrade Soil Equilibrium.—As shown Fig. the pressure applied 
the wheels and aggregate layer the portion the geotextile given 


Eq. 14. The pressure applied the subgrade soil the portion 
the geotextile is: 


which: reduction pressure resulting from the use geotextile 
(as explained previously) Pascals; will evaluated the next sections. 

Because the subgrade soil confinement provided the geotextile helps keep 
the deflection small for all applied pressures less than the ultimate bearing 
capacity, the pressure can large the ultimate bearing capacity 
the subgrade soil, given Eq. 16, which replaced yh: 


comparison between and outlines the benefit resulting from 
subgrade soil confinement geotextile; the term Eq. referring 


sets of dual wheels 


aggregate 


eer 


geotextile 
c 


subgrade soil 


FIG. 10.—Kinematics Unpaved Road with Geotextile 


elastic behavior replaced Eq. the term 2)c, referring 
plastic behavior. 


previously stated, function the tension the geotextile, which 
depends its elongation, the elongation itself depending its shape. Therefore, 
prior the determination p,, the shape and the resulting elongation the 
geotextile must assessed; this done the following two sections. 

Shape Deformed wavy shape the deformed geotextile 
results from the incompressibility the subgrade soil. Therefore, the volume 
subgrade soil displaced downwards settlement must equal the volume 
soil displaced upwards heave. The shape the deformed geotextile 
assumed consist portions parabolas connected and points located 
the initial plane the geotextile (Fig. 11). also assumed that the thickness 
aggregate layer not significantly affected subgrade soil deflection. This 
assumption based observations from full scale tests (4). Therefore: (1) 
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the total deflection the geotextile equal rut depth; and (2) comparison 


Two cases must considered, depending upon which parabola the wider. 
The first case arises when the central parabola the wider (a’ thus, 
according Eqs. and 31: tan 2B). One half the soil displaced 
downwards settlement contributes heave between the wheels, while the 
other half contributes heave beyond the wheels [Fig. 12(a)]. The relationship 
between settlement, and rut depth, which characterizes the shape the 
geotextile determined writing that areas and are equal [Fig. 12(a)]: 


Therefore: 


The second case arises when the parabolas under the wheels are the wider 
a’). Less than one half the soil displaced downwards settlement 


of geotextile 


FIG. 11.—Shape Deformed Geotextile 


contributes heave between the wheels [Fig. 12(b)]. The width area 
assumed equal half the width (P’) and the relationship between 
and determined writing that areas and are equal [Fig. 
After calculations: 


2 


2ra 


Note that Eqs. and lead the same value, r/2, when a’. 

both Eqs. and 34, and could replaced their expressions 
from and order obtain the relationship between and 

Elongation Geotextile.—The reduction pressure, p,, due the tension 
geotextile parabola (P) (Fig. 12). Since tension function elongation, 
necessary determine first the elongation, the geotextile (P). 
Again, the same two cases must considered. 

elongation, and the tension the geotextile tend larger (P) than 


(30) 
2as 
ra’ 
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(P’). result, the geotextile (P’) pulled towards (P). most cases, 
the overburden over the portion (P’) the geotextile too small generate 
sufficient friction prevent movement, and the geotextile actually moves from 
(P’) (P). Therefore, reasonable assume that elongation uniform 
along the entire length the geotextile: 


which: elongation the geotextile; half length the chord 
subtended respectively (P) and (P’) (m); and half length (P) 
and (P’), respectively (in meters). 


initial location 
of geotextile 


geotextile 


initial location 
of geotextile 


FIG. 12.—Determination Settlement Versus Rut Depth Relationship 


The elongation can calculated using Eq. and the following relationship 
between arc parabola and subtended chord: 


2 , 


elongation (consequently the tension) the geotextile larger (P’) than 
(P). However, the geotextile does not move because the high normal stress 
transmitted (P) the wheel load generates important friction. So, 
this case, the elongation different (P) and (P’). The elongation the 
geotextile (P) given directly Eq. because: 


a+a 
area A 
a 
a 
a 
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Determination Reduction Pressure Subgrade Soil.—In fact, 
uniform pressure applied (Fig. 11) and equivalent the vertical 
projection the tension the geotextile points and 


Determination Thickness Aggregate Layer.—The relationship between 


cohesion subgrade soil, c,, and thickness aggregate layer, stems from 
Egs. 14, 29, and 42: 


are given respectively Eqs. and (on-highway trucks) and and 
(off-highway trucks); and and are given respectively Eqs. and 
all cases. 

Fig. presents values obtained from Eq. the case the standard 
axle load and 0.3 rut depth. Similar charts have been established for several 
values axle load and rut depth Design engineers should bear mind 
that the use these charts restricted very light traffic 
passages, approx) since has been established using quasi-static analysis. 
fact, better consider these charts step towards the more general 
method presented hereafter. 


Principle Method.—To this point, unpaved roads without and with geotextile 
have been considered. For unpaved roads without geotextile, two approaches 
have been presented: (1) quasi-static analysis leading thickness 
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According property parabolas: 
a 
From Fig. 
Combining 39, 40, and 41: 
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aggregate layer; and (2) empirical method taking traffic into account, giving 
thickness aggregate layer. For unpaved roads with geotextile, only 
quasi-static analysis has been presented, leading thickness aggregate 
layer. 


--- tire inflation pressure * 480 kPa 
—— tire inflation pressure os 620 kPa 


450 kN/m 

400 kN/m 

K = 300 kN/m geotextile 
200 kN/m modulus 
100 


K = 10 kKN/m 


(kPa) 


CBR 


FIG. 13.—Aggregate Thickness Versus Subgrade Soil Cohesion (Quasi-Static Analy- 
sis for Case with Geotextile) 


The thickness aggregate layer when traffic taken into account remains 
determined the case unpaved roads with geotextile. The following 
procedure proposed: 


The reduction aggregate thickness, Ah, resulting from the use 
geotextile, deduced from the quasi-static analyses by: 


The thickness, h’, aggregate layer the case geotextile-reinforced, 
unpaved road when traffic taken into account, determined by: 


This procedure assumes that the value does not depend the traffic. 
Extensive full scale tests should performed order verify this assumption. 
Practical Charts.—The establishment charts giving the thickness, h’, 
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function all the parameters would too cumbersome. way minimizing 
the number charts take advantage the fact that does not depend 
the geotextile while does not depend the traffic. Therefore, two sets 
curves are presented each chart, one giving function the traffic, 
the other giving function the geotextile modulus. Then, 
the user the chart make subtraction, according Eq. 45, determine 
the aggregate thickness, h’, taking traffic and geotextile modulus into account. 
example such chart given Fig. 14. The values are extracted 
from Fig. and the values are deduced from Figs. and using Eq. 
44. Fig. also presents set curves giving the geotextile elongation versus 
subgrade soil cohesion c,. This set curves was obtained combining 


=P,= 80 KN 
=0.3m 
P, = 480 kPa 


Sn for 


K 450 
K = 400 
K © 300 KN/m 
K = 200 kN/m 
K = 100 KN/m 
10 kN/m 


elongation 


of geotextile ——+ 


number of 


passages 


FIG. 14.—Aggregate Thickness (Case without Geotextile when Traffic Taken 
into Account) and Possible Reduction Aggregate Thickness Resulting from 
use Geotextile Versus Subgrade Soil Cohesion (Chart Related On-Highway 
Truck with Standard Axle Load) 


30-38 and Eq. (for on-highway trucks) Eq. (for off-highway trucks). 
allows the user the chart check that, the considered case, the geotextile 
not subjected excessive elongation. 

Two features the chart are noteworthy: (1) according Eq. 44, can 
never higher than this the reason why all the versus curves 
end the versus extracted from Fig. and (2) aggregate needed 
top the geotextile when the versus curve above the versus 
curve. 

Analysis Results.—Fig. only one among the charts developed for 
various axle loads and rut depths. Considering the whole set charts leads 
the following comments: (1) the considered range tire inflation pressures, 


(m) 
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the tire inflation pressure has influence and almost none Ah; (2) 
for typical geotextiles, the reduction aggregate thickness resulting from the 
use geotextile generally ranges between 20% and 60% for subgrade soil 
CBR one and number passages between 1,000 and 10,000; however, 
the use geotextile with high modulus makes possible build access 
road designed for very small number passages, with aggregate all 
(this may useful for military applications for soil exploration); (3) 
increase the allowable rut depth causes decrease (thickness aggregate 
layer without geotextile) and increase (reduction thickness resulting 
from the use geotextile); and (4) the higher the geotextile tensile modulus, 
the thinner the required thickness aggregate layer. 

Practical Example.—The example presented Appendix shows how given 
traffic can replaced 340 passages single axle load about kN. 


TABLE 1.—Full Scale Test Data Versus Calculated Aggregate Thickness 


Aggre- 
gate Thick- 
ness (for 
Rut Depth 


FULL SCALE TEST DATA (FOR RUT DEPTH 0.28 
Geotextile 


from 
CBR, and 


kilonewtons Passages 


per meter CBR meters meters 
(3) (4) 


Without 
geotextile 

Bidim U44/C38 

Reeves T16 


values because number passages larger than 10,000. 


What the required thickness aggregate layer, the CBR the subgrade 
soil one, the modulus the geotextile and the allowable rut 
depth 0.3 

Fig. 14, related axle load and 0.3 rut depth: (1) Point 
corresponding CBR and 340, shows that the required thickness 
aggregate when there geotextile would 0.35 and (2) point 
corresponding CBR and kN/m, shows that the reduction 
thickness can 0.15 Consequently, the recommended thickness 
aggregate 0.20 

Also, point Fig. shows that the elongation the geotextile approx 
10%. should checked that the elongation failure the geotextile 
larger than this value. 

Comparison with Full Scale Tests.—The United States Army Corps Engineers 


(1) 
200 0.36 0.34 
2,500 0.36 0.38 
200 37,000 0.36 
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has conducted full scale tests unpaved roads with and without geotextile 
(3,4). The results these tests were expressed number passages 
equivalent standard axle (P, kN). So, comparison easy with the 
theoretical results presented Fig. 14. Test results corresponding rut depth 
0.28 are given Table Also tabulated the aggregate thickness calculated 
from test data using Fig. 14. 

Comparison calculated and actual thicknesses shows good agreement 
when traffic light. The theoretical results appear conservative when traffic 
heavy. This confirms that the use Fig. and similar charts not 
recommended for numbers passages larger than 10,000. 


Summary Main Features this Study.—The development method 
for the design geotextile-reinforced, unpaved roads involves three steps: 


empirical formula deduced from full scale test results related 
unpaved roads without geotextile; gives the thickness aggregate layer 
function traffic and soil properties. 

theoretical analysis compares unpaved roads with and without geotextile; 
determines the reduction aggregate thickness resulting from the use 
geotextile; this theoretical analysis does not take traffic into account. 

the empirical formula and the theoretical analysis, charts 
are established; these charts allow the determination aggregate thickness 
for unpaved road with geotextile, when traffic taken into account. 


The method presented this paper devoted the reinforcement role 
geotextiles and therefore, does not consider other beneficial effects 
geotextiles, such separation, filtration, and drainage. 

Validity Assumptions.—Many assumptions were made this study. All 
are precisely defined the text. However, short review necessary remind 
the potential users the charts their limitations: 


Assumptions related bearing capacity calculations are classical soil 
mechanics (bidimensional calculations, use the theories elasticity and 
plasticity); therefore, from the viewpoint bearing capacity, there are more 
limitations than, say, classical methods for foundation design. 

Assumptions related the role the geotextile are based upon limited 
experimental evidence; these assumptions may improved the future; but 
all assumptions, made this study, concerning the mechanical roles soil 
and geotextile are consistent with each other. 

Assumptions related traffic are the most imprecise; the fact that Eq. 
(established for paved roads) used for unpaved roads can questioned, 
but, present, there similar equation for unpaved roads; safe, 


recommended not use Fig. and similar charts for numbers passages 
larger than 10,000. 


Practical Consequences the Results.—The present charts are intended 
assist engineers designing geotextile-reinforced, unpaved roads and particu- 
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lar, evaluating the reduction aggregate thickness resulting from the use 
geotextile. This reduction appears the order magnitude 
20%-60% some typical cases considered this paper. Such significant 
reduction can have important economic consequences, that justify the use 
the charts. 

These charts, related frictionless subgrade soils, have already proved 
useful several projects. hoped that further research will lead 
similar charts for other types subgrade soils. 
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Axte 


traffic defined the load, P,, and the number passages, N,, 
each axle. order use the charts presented this paper, traffic must 
expressed the simplest manner which consists the number passages, 
equivalent single axle load, This can achieved using 


Eq. and assuming that all axles have the same width and are driven 


the same ruts (in other words, the trucks are assumed not wander over 
the road). 


General Case.—There are many possibilities define the equivalent single 
axle corresponding given traffic. Two these possibilities seem the most 
adequate. 


first possibility consider that the load the equivalent axle equal 
the largest 


(P,) max 
Hence, from Eq. 25: 


3.95 


second possibility consider that the number passages the equivalent 
single axle equal the largest 


(N)max 
Thus, from Eq. 25: 


These two approaches are logical because damage the road can caused 
either excessive loading excessive number passages. 
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Then, charts are used with both and and the more 
conservative result selected the design thickness aggregate layer. 

Case Identical Trucks with Rear Single Axle.—Consider passages 
identical trucks with rear single axle. Because the rear axle load larger 


than the front axle load, Eqs. and become: 


tf 


which: front axle load and rear axle load (N). Usually, 
least two times larger than Consequently: 


eq 

So, practically only the rear axle needs considered when the traffic 
made identical trucks with rear single axle. 

Case Identical Trucks with Rear Tandem Axle.—According the previous 
comment, the front axle load, smaller than the rear axle load, can neglected. 
Therefore, only the two axles the tandem need considered. Each 
these two axles passing times the road: 


The equivalent single axle load should Q/2 (in which tandem axle 
load). However, recommended take: 


This increase 20% the actual equivalent single axle load intended 
take into account the additional stresses induced the subgrade soil because 
the interaction between the two axles the tandem. 

Example Determination Equivalent Single Axle.—Consider passages 
trucks with rear tandem axle load 160,000 and 340 passages trucks 
with rear single axle load 50,000 What the equivalent single axle 
used the design the road? 

First, according Eqs. and 55, the tandem axle equivalent single 
axle defined by: 


Then, there are two alternatives determine the single axle equivalent 
the total traffic. The first alternative is, according Eqs. and 47: 


50,000 
96,000 


The second alternative is, according Eqs. and 49: 


140 1/3.95 
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Then, the aggregate thickness calculated using the charts for: (1) 166 passages 
single axle 96,000 and (2) 340 passages single axle 80,047 
The more conservative result must selected. 
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The following symbols are used this paper: 


area shown Fig. 12, square meters; 
area shown Fig. 12, square meters; 
area shown Fig. 12, square meters; 
contact area tire, square meters; 
half length chord subtended portion parabola (P), 
shown Fig. 11, meters; 
half length chord subtended portion parabola (P’), 
shown Fig. 11, meters; 
width equivalent contact area, meters; 
half length portion parabola (P), meters; 
half length portion parabola (P’), meters; 
California Bearing Ratio (dimensionless); 
undrained cohesion subgrade soil (in pascals kilopascals); 
distance between midpoints two sets wheels, meters; 
thickness subgrade soil, meters; 
minimum value ensure validity analysis, meters; 
depth plastic zone, meters; 
aggregate thickness (quasi-static analysis case with geotextile), 
meters; 
aggregate thickness (case with geotextile, when traffic taken 
into account), meters; 
aggregate thickness (quasi-static analysis case without geo- 
textile), meters; 
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aggregate thickness (case without geotextile, when traffic taken 
into account), meters; 
reduction aggregate thickness from use geotextile, meters; 
secant modulus geotextile (in Newtons per meter and kilonew- 
ton per meter); 
number passages axle load 
equivalent number passages (Appendix 
number passages standard axle load 
axle load (in Newtons and kilonewtons); 
equivalent single axle load (Appendix I), pascals kilopascals; 
front axle load (in Newtons and kilonewtons); 
rear single axle load (in Newtons and kilonewtons); 
standard axle load (in Newtons and kilonewtons); 
pressure base aggregate layer case with geotextile 
pascals kilopascals; 
pressure base aggregate layer case without geotextile 
pascals kilopascals; 
pressure applied subgrade soil portion geotextile 
pascals kilopascals; 
tire inflation pressure pascals kilopascals; 
equivalent contact pressure pascals kilopascals; 
reduction pressure resulting from use geotextile pascals 
kilopascals; 
tandem axle load (in Newtons and kilonewtons); 
cone resistance pascals kilopascals; 
elastic bearing capacity soil pascals kilopascals; 
lateral surcharge pascals kilopascals; 
ultimate bearing capacity soil pascals kilopascals; 
rut depth, meters; 
settlement under wheels, meters; 
tension portion geotextile, shown Fig. (in 
Newtons per meter kilonewtons per meter); 
angle load distribution pyramide case with geotextile, 
degrees; 
angle load distribution pyramide case without geotextile, 
degrees; 
angle between tension and vertical line, shown Fig. 11, 
degrees; 
weight per unit volume aggregate (in Newtons per cubic meters 
kilonewtons per cubic meters); 
elongation geotextile parabola (P) (dimensionless, 
percentage); and 

elongation geotextile failure (dimensionless, percentage). 


(Note: this paper, elongation synonymous strain) 
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